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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non— SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply 

cubic  feet 
feet 


_ Bx _ 

0.02831685 

0.3048 


To  Obtain 
cubic  metres 
metres 
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USER'S  GUIDE  FOR  A  THREE-m MEWS TONAL  NUMERICAL 


HYDRODYNAMIC.  SALINITY.  AND  TEMPERATURE  MODET. 

OF  CHESAPEAKE  BAY 

PART  I :  INTRODUCTION 

1.  Chesapeake  Bay  is  one  of  the  largest  estuaries  in  the  world.  As 
shown  in  Figure  1,  the  main  bay  extends  approximately  300  km  north  from  the 
ocean  entrance  to  the  Susquehanna  River.  The  average  depth  of  the  bay  is 
about  8  m,  although  a  natural  channel  with  depths  greater  than  15  m  traverses 
the  bay  for  more  than  60  percent  of  its  length.  The  bay  is  irregular  in 
shape,  varying  in  width  from  6.4  km  to  48.3  km  in  the  middle  bay  off  the 
Potomac  River.  The  bay  is  sufficiently  long  to  accommodate  one  complete  tidal 
wave  at  all  times . 

2.  Chesapeake  Bay  is  one  of  the  Nation's  most  valuable  natural 
resources.  It  supports  Important  commercial  and  recreational  fisheries, 
transportation,  industry,  recreation,  and  tourism,  and  provides  irreplaceable 
habitat  for  living  marine  resources  and  wildlife.  However,  the  estuary  has 
been  subjected  to  increasing  environmental  stress  in  recent  decades,  and  the 
productivity  and  beauty  of  Chesapeake  Bay  have  significantly  declined. 

3.  In  1983,  the  US  Environmental  Protection  Agency  (USEPA)  identified 
major  contributing  factors  to  the  bay's  decline  as  inputs  of  nutrients  and 
toxicants  from  point  and  nonpoint  sources,  changes  in  land  use  within  the 
basin  with  resulting  modification  of  the  watershed  and  resource  habitat,  and 
concurrent  impacts  of  natural  events  such  as  floods  and  droughts  (USEPA  1983a, 
1983b).  Because  population  within  the  bay  drainage  basin  is  still  increasing, 
and  development  pressures  will  continue  into  the  forseeable  future,  it  is 
necessary  that  strategies  are  developed  to  reverse  the  present  bay  decline  and 
to  accommodate  future  growth  in  an  environmentally  sound  manner. 

4.  The  Chesapeake  Bay  Program  (CBP),  established  in  1983,  provides  a 
management  structure  through  which  the  activities  of  state  and  Federal  agen¬ 
cies,  as  well  as  those  of  private  citizens,  can  be  coordinated  toward  the  goal 
of  bay  restoration.  A  number  of  tools  have  been  developed  by  the  CBP  which 
will  assist  in  the  planning,  implementation,  and  evaluation  of  these  strate¬ 
gies.  Among  these  are  a  three-dimensional  (3D),  time-varying  water  quality 
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model  for  the  bay  and  its  tributaries  to  provide  a  detailed  assessment  of  the 
system's  response  to  nutrient  inputs  and  other  parameters  varied  realistically 
over  time  and  space.  The  3D  water  quality  model  represents  the  means  through 
which  proposed  management  actions  can  be  tested  before  implementation,  thereby 
allowing  more  cost-effective  selection  of  appropriate  strategies.  However,  it 
has  been  acknowledged  that  water  quality  impacts  in  Chesapeake  Bay  cannot  be 
successfully  assessed  without  an  accurate  description  of  the  hydrodynamic 
processes.  Thus,  a  3D  time-varying  hydrodynamic  model  of  Chesapeake  Bay  has 
been  developed  at  the  US  Army  Engineer  Waterways  Experiment  Station  (WES) . 

The  model  is  called  CH3D— WES  (Curvilinear  Hydrodynamics  in  Three  Dimensions) . 

5.  The  numerical  grid  employed  in  the  Chesapeake  Bay  3D  hydrodynamic 
model  is  shown  in  Figure  2.  There  are  734  active  horizontal  cells  and  a  maxi¬ 
mum  of  15  vertical  layers,  resulting  in  3992  computational  cells.  To  capture 
the  important  features  of  the  hydrodynamic  processes  and  bathymetry  in  the 
bay,  grid  resolution  is  1.52  m  vertical  and  approximately  10  km  longitudinal 
and  3  km  lateral.  Major  tributaries  are  modeled  fully  3D  in  the  lower  reach 
and  as  constant-width  two-dimensional  (2D)  in  the  upper  reach. 

6.  This  user's  guide  presents  a  detailed  discussion  of  theoretical 
aspects  of  the  3D  model.  The  governing  equations  solved  in  the  3D  model  are 
presented  in  both  Cartesian  and  curvilinear  or  boundary-fitted  form  for  both 
the  external  and  internal  modes.  The  external  mode  consists  of  vertically 
averaged  equations,  whereas  the  internal  mode  contains  the  full  3D  equations. 
Required  boundary  conditions,  along  with  a  discussion  of  the  algorithm 
employed  to  numerically  solve  the  equations,  are  given.  In  addition,  the 
structure  of  the  computer  model  is  discussed.  Detailed  descriptions  of  the 
function  of  each  subroutine  are  provided  along  with  a  discussion  of  the  tiow 
of  c-imputations  in  the  model.  The  type  of  input  data  required  to  operate  the 
model  and  the  setup  of  those  data  are  also  presented. 

7.  Basic  input  required  are  time-varying  water-surface  elevations  at 
the  ocean  entrance  and  freshwater  inflows  at  the  head  of  all  tributaries.  In 
addition,  time-varying  salinity  and  temperature  data  must  be  prescribed  at  all 
inflow  boundaries.  Also,  time-varying  wind  and  surface  heat  exchange  data 
must  be  prescribed  at  one  or  more  locations.  These  are  then  distributed  in 
some  manner  over  the  entire  grid.  All  input  data,  including  initial  condi¬ 
tions,  bathymetry,  boundary,  and  computational  control  data  are  input  from 
fixed  files.  An  application  using  data  from  September  1983  serves  to 
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demonstrate  the  setting  up  of  input  files.  A  detailed  list  of  the  input  data 
requirements  and  the  proper  format  are  also  given.  Standard  output  are  water- 
surface  elevations,  velocity,  salinity,  and  temperature  in  each  cell  of  the 
grid  at  some  specified  time  interval. 


9 


PART  II:  THEORETICAL  BASIS  OF  CHESAPEAKE  BAY 
3D  HYDRODYNAMIC  MODEL 


8.  As  noted,  the  particular  hydrodynamic  computer  model  employed  is 
called  CH3D-WES.  The  basic  model  (CH3D)  was  developed  by  Sheng  (1986)  for  WES 
but  has  been  extensively  modified  in  its  application  to  Chesapeake  Bay.  These 
modifications  have  consisted  of  different  basic  formulations  as  well  as  sub¬ 
stantial  recoding  for  more  efficient  computing.  As  its  name  implies,  CH3D-WES 
makes  hydrodynamic  computations  on  a  curvilinear  or  boundary-fitted  planform 
grid.  Physical  processes  impacting  bay— wide  circulation  and  vertical  mixing 
that  are  modeled  include  tides,  wind,  density  effects  (salinity  and  tempera¬ 
ture),  freshwater  inflows,  turbulence,  and  the  effect  of  the  earth's  rotation. 
Adequately  representing  the  vertical  turbulence  is  crucial  to  a  successful 
simulation  of  stratification/destratification.  A  second— order  turbulence 
model  based  upon  the  assumption  of  local  equilibrium  of  turbulence  is 
employed.  The  boundary— fitted  coordinates  feature  of  the  model  provides 
enhancement  to  fit  the  deep  navigation  channel  and  irregular  shoreline  config¬ 
uration  of  the  bay  and  permits  adoption  of  an  accurate  and  economical  grid 
schematization.  The  solution  algorithm  employs  an  external  mode  consisting  of 
vertically  averaged  equations  to  provide  the  solution  for  the  free  surface  to 
the  internal  mode  consisting  of  the  full  3D  equations. 


Basic  Equations 


9.  The  basic  equations  for  an  incompressible  fluid  in  a  right-handed 
Cartesian  coordinate  system  (x,y,z)  are: 


3u  3v  9w  p, 

Sy  -Ji 


(1) 


9u  du^  9uv  9uw  _  1  3p  d  Ia  3u 

■?t  Sy  ^ 
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(6) 


P  =  P  (T.S) 


(7) 


where 

(u,v,w)  «  velocities  in  (x,y,z)  directions 
t  -  time 

f  -  Coriolis  parameter  defined  as  20  sin  ^ 
where 

O  -  rotational  speed  of  the  earth 
^  -  latitude 

p  -  density 

p  »  pressure 

A„,K„  -  horizontal  turbulent  eddy  coefficients 
Ay.Ky  -  vertical  turbulent  eddy  coefficients 

g  -  gravitational  acceleration 

T  -  temperature 

S  -  salinity 

Equation  4  implies  that  vertical  accelerations  are  negligible  and  thus  the 
pressure  is  hydrostatic. 
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10.  Various  forms  of  the  equation  of  state  can  be  used  for  Equation  7. 
In  the  present  model,  the  equation  given  below  is  used; 


p  =  P/(a  +  0.698P) 


(8) 


where 


P  =  5890  +  38T  -  0.375T2  +  3S 


a  =  1779.5  +  11.25T  -  0.07451^  -  (3.8  ♦  0.01T)S 


and  T  is  in  degrees  Celsius,  S  is  in  parts  per  thousand,  and  p  is  in 
grams  per  cubic  centimeter. 

Non-Dimensionalizatlon  of  Equations 

11.  Working  with  the  dimensionless  form  of  the  governing  equations 
makes  it  easier  to  compare  the  relative  magnitude  of  various  terms  in  the 
equations.  Therefore,  the  following  dimensionless  variables  are  used; 


(u*,  v*.  w*)  X  (u.  V,  wX,/Z,)/U, 
(x‘,  y-,  z*)  X  (X.  y.  zX,/Z.)/X, 
’■;)  =  (»■;.  ’■p/PofZ.U. 

t*  X  tf 

r*  X  gf/fu,x,  X  f/s. 

P*  =  (p  -  -  Po) 


T*  X  (T  -  T„)/(T,  -  T„) 
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K  =  A^/Avr 
K  =  K„/K„r 
k;  =  K/Kr 


where 

(r*,  r"  )  -  wind  stress  in  (x,y)  directions 

X  y 

f  “  water  surface  elevation 

which  then  yields  the  following  dimensionless  parameters  in  the  governing 
equations : 


Vertical  Ekman  Number: 

K  - 

=  A,,/fZ,2 

Lateral  Ekman  Number: 

Eh  = 

AHr/fXr" 

Vertical  Prandtl  (Schmidt)  Number: 

=  Ki/^z 

Lateral  Prandtl  (Schmidt)  Number: 

~  Anr/^r 

Froude  Number: 

Fr  = 

=  U,/(gZ,)i/2 

Rossby  Number: 

Bo  = 

=  VfXr 

Densimetric  Froude  Number: 

Fr„ 

=  F,  /yfT 

where 

c  = 

<Pr  -  Po)/Po 

Sr  .  Tj  ,  Uj  ,  Pr  ,  Xj  ,  ,  A„t  ,  ,  Khj  ,  and  are  arbitrary 

reference  values  of  the  salinity,  temperature,  velocity,  density,  etc. 
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External-Internal  Modes 


12.  The  basic  equations  (Equations  1-4)  can  be  integrated  over  the 
depth  to  yield  a  set  of  vertically  integrated  equations  for  the  water  surface , 
f  ,  and  unit  flow  rates  U  and  V  in  the  x  and  y  directions.  Using  the 
dimensionless  variables  (asterisks  have  been  dropped)  and  the  parameters  pre¬ 
viously  defined,  the  vertically  integrated  equations  constituting  the  external 
mode  are : 


=  -  H  4£  .  r.,  -  r,,  .  V 

Tt  3x 


-  R. 


^  ^  Is]  *  ^  l|]] 


_  Rd  dp 
Fr^ 


av  _  H  +  r  -  r  -  U 

TE  Ty 


^  ^  Is]  *  l|]] 


where 


_  Rd  H2  ap 


fi  =  gZr/f2X,^  -  (R„/F,)2 


(11) 


(12) 


(13) 
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As  will  be  discussed  later,  the  major  purpose  of  the  external  mode  is  to  pro¬ 
vide  the  updated  water-surface  field. 

13.  The  internal  mode  equations  from  which  the  3D  velocity,  salinity, 
and  temperature  fields  are  computed  are: 


K  ^  L 

dhul 

FSi 

^  -Ji 

"FT] 

dhuu 

dhuv 

dhuw 

dx 

*  ST  * 

dz 

s 

-  ^  (/:  ^  H 

P  f3hvu  ^  dhw  ^  9hvw  I 

-3y-  -FT  J 


(14) 


(15) 


duh  ^  dvh 
5x  5y 

L  J 


(16) 


dhT  _  ^  5  fi;  3t1  _  p  fdhuT  dhvT  dhwT 

■  T?;  -Ji  •jiJ  ■  "  ST  " 


E,, 

TF 


d 

’  dhl' 

d 

V  ^hll 

~F^ 

*  ^ 

(17) 


(18) 


ahs  _  E*,  a  as]  _  p  fahas  ^  ahvs  ^  ahws] 
at  Pr^  "Jz  dx  ay  dz  J 


In  these  equations  h  is  the  thickness  of  an  internal  layer  and  k+1/2  and 
k-1/2  represent  the  top  and  bottom,  respectively,  of  the  k^*'  vertical 
layer . 


Boundary-Fitted  Equations 

14.  To  better  resolve  complex  geometries  in  the  horizontal  directions, 
the  Chesapeake  Bay  3D  model  makes  computations  on  the  boundary-fitted  or  gen¬ 
eralized  curvilinear  planform  grid  shown  in  Figure  2.  This  necessitates  the 
transformation  of  the  governing  equations  into  boundary-fitted  coordinates 
(i,Ti).  If  only  the  (x,y)  coordinates  are  transformed,  a  system  of  equations 
similar  to  those  solved  by  Johnson  (1980)  for  vertically  averaged  flow  fields 
is  obtained.  However,  in  the  Chesapeake  Bay  3D  model  not  only  are  the  (x,y) 
coordinates  transformed  into  the  (^,i?)  curvilinear  system  but  also  the 
velocity  is  transformed  such  that  its  components  are  perpendicular  to  the 
(^,17)  coordinate  lines.  This  is  accomplished  by  employing  the  definitions 
below  for  the  components  of  the  Cartesian  velocity  (u,v)  in  terms  of 
contravariant  components  u  and  v 


U  =  XjU  +  x„v 

V  =  y^u  +  y^-v  (19) 


along  with  the  following  expressions  for  replacing  Cartesian  derivatives 

t.  ■  4  [  (*y-)e  -  1 

f,  .  1  [-  (fx,),  .  (fx<),  ] 
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where  J  is  the  Jacobian  of  the  transformation  defined  as 


J  -  ^yt 

With  the  governing  equations  written  in  terms  of  the  contravariant  components 
of  the  velocity,  boundary  conditions  can  be  prescribed  on  a  boundary-fitted 
grid  in  the  same  manner  as  on  a  Cartesian  grid  since  u  and  v  are  perpen¬ 
dicular  to  the  curvilinear  cell  faces  (e.g.,  at  a  land  boundary,  either  u 
or  V  is  set  to  zero) . 

15.  Initially  the  vertical  dimension  was  handled  through  the  use  of 
what  is  commonly  called  a  sigma-stretched  grid.  Such  a  grid  is  illustrated  in 
Figure  3.  However,  with  the  grid  resolution  shown  in  Figure  2  it  was  observed 
that  stratification  in  the  deep  channels  could  not  be  maintained  during  long 
term  simulations  on  Chesapeake  Bay.  With  a  sigma-stretched  grid,  the  bottom 
layer  in  one  column  communicates  with  the  bottom  layer  in  an  adjacent  column. 
Thus,  if  depth  changes  are  rather  coarsely  resolved,  channel  stratification 
cannot  be  maintained.  As  a  result,  the  governing  equations  presented  below 
for  solution  on  the  Cartesian  or  z-plane  in  the  vertical  direction  are  the 
ones  constituting  the  internal  mode  in  the  Chesapeake  Bay  model. 


cr  = 


ili 


Figure  3.  Sigma-stretched  grid 
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16.  With  both  the  Cartesian  coordinates  and  the  Cartesian  velocity 
transformed,  the  following  boundary-fitted  equations  for  u  ,  v  ,  w  ,  S 
and  T  to  be  solved  in  each  vertical  layer  are  obtained. 


■5? 


(Jy^h 


lUU 


+  Jy„huv  )  +  ^  (Jy^huv  +  Jy,hv^|  -  ^2  1.^  (Jx^huu  +  Jx„hu^ 


a 


(  JXfhuv  +  Jx,hvv  )j  -  R®  [(wu  )top  -  (wu  )i,^J 

C  au1  _  C  au]  _  f  p 


f  p22  dp 

■?? 


dz  +  Horizontal  Diffusion 


^12  Spl 

■  ■3^J 

+  Jy„hu^  +  ^  (Jy^huv  +  Jy^hv^j  +  ^  (Jx^huu  +  Jx^hH^ 


a  —  — 

+  .g—  (  JXjhuv  +  Jx,,hvv ) 


R®  [(wv)t®p  -  (wvLt 


+  Horizontal  Diffusion 


(21) 


(22) 


«  =  -  1  faJuh  ^  aJvh] 

Wtop  Wb®t  3 


(23) 
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ahs  asl  _  asl  _  fahJus  ^  ahjvs’) 

-3F  =  -51^  ~  T  [-5?-  -5^J 


-  Ro  [(wS)top  -  (WS)^„J  ^  Horizontal  Diffusion 


ahT 

-  ^ 

<,3X1  _  [ 

,  ail 

R. 

fahJuT  ^ 

ahjvrl 

■JF 

■  ^  [ 

'”■5^1 

h?- 

-  Rp  [(wT)t„p  -  (wT)j^t]  *  Horizontal  Diffusion 

where 

Gii  =  x|  +  yf 
G22  =  ^  y" 

Gi2  =  G21  =  +  y^y, 

17.  Similarly,  the  transformed  external  mode  equations  become: 


(24) 


(25) 


(26) 


(27) 


au  H  L  ar  p  ar 

■5t  “  ■  P2  -3^  -  Gi2  -3^ 


1  (G12U  -  G22^  >  ^  (jy^UU  .  Jy,uv)  .  ^  (jy,uv  .  Jy^Vv) 


Rt,y„ 


■5? 


^  ’•.f  -  ’■bf 


(  JXfUU  +  Jx„U^  + 

a 

(jXjUV 

Ro  h2 

Ig  - 

dp\ 

'^12 
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+  Horizontal  Diffusion 


(28) 


av 

■Ft  “ 


T 


+  Horizontal  Diffusion 


(29) 


18.  Equations  27-29  are  solved  first  to  yield  the  water-surface  eleva¬ 
tions  which  are  then  used  to  evaluate  the  water-surface  slope  terms  in  the 
internal  mode  equations.  The  horizontal  diffusion  terms  are  quite  lengthy  and 
thus  are  presented  separately  in  Appendix  A  for  the  internal  mode  equations. 
Similar  expressions  for  the  vertically  averaged  equations  can  be  inferred  from 
those  for  the  internal  mode . 


Boundary  Conditions 

19.  The  boundary  conditions  at  the  free  surface  are 


3T 

Ti 


Pr 


K  (T  -  TJ 


as 

■F5 


=  0 


(30) 
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whereas  the  boundary  conditions  at  the  bottom  are 


-i) 


ax  as  _  0 

Si  '  SI 


(31) 


where 

C  -  surface  drag  coefficient 
W  -  wind  speed 

K  -  surface  heat  exchange  coefficient 
T,  -  equilibrium  temperature 
C<j  -  bottom  friction  coefficient 

Ui  ,  Vj  -  values  of  the  horizontal  velocity  components  next  to  the 
bottom. 

With  Zj  equal  to  one-half  the  bottom  layer  thickness,  Cj  is  given  by 


Cj  *  k2  [in  (zj/z^)] 


(32) 


where 

k  -  von  Karman  constant 
Zo  -  bottom  roughness  height 

Manning's  formulation  is  employed  for  the  bottom  friction  in  the  external  mode 
equations  if  the  model  is  used  purely  to  compute  vertically  averaged  flow 
fields.  As  presented  by  Garratt  (1977),  the  surface  drag  coefficient  is 
computed  from 


C  -  (0.75  ♦  0.067  W)  X  lO'® 


(33) 
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with  the  maximum  allowable  value  being  0.003.  The  surface  heat  exchange  coef¬ 
ficient,  K  ,  and  the  equilibrium  temperature,  T„  ,  are  computed  from 
meteorological  data  (wind  speed,  cloud  cover,  wet  and  dry  bulb  air  tempera¬ 
tures,  and  relative  humidity)  as  discussed  by  Edinger,  Brady,  and  Geyer 
(1974) . 

20.  Along  the  shoreline  where  river  inflow  occurs,  the  freshwater  in¬ 
flow  and  its  temperature  are  prescribed  and  the  salinity  is  assumed  to  be 
zero.  At  an  ocean  boundary,  the  water-surface  elevation  is  prescribed  along 
with  time-varying  vertical  distributions  of  salinity  and  temperature.  During 
flood,  the  specified  values  of  salinity  and  temperature  are  employed,  whereas 
during  ebb,  interior  values  are  advected  out  of  the  grid.  Along  a  solid 
boundary,  the  normal  component  of  the  velocity  and  the  viscosity  and 
diffusivity  are  set  to  zero. 


Initial  Conditions 

21.  Uhen  initiating  a  run  of  the  Chesapeake  Bay  model,  the  values  of 
f  ,  u  ,  V  ,  w  ,  U  and  V  are  all  set  to  zero.  Values  of  the  salinity 
and  temperature  are  read  from  input  files.  These  initial  fields  are  generated 
from  known  data  at  a  limited  number  of  locations.  Once  the  values  in  individ¬ 
ual  cells  are  determined  by  interpolating  from  the  field  data,  the  resulting 
3D  field  is  smoothed  several  times.  Generally,  the  salinity  and  temperature 
fields  are  frozen  for  the  first  few  days  of  a  simulation. 

Numerical  Grid 

22.  A  staggered  grid  is  used  in  both  the  horizontal  and  vertical  direc¬ 
tions  of  the  computational  domain  (Figure  4).  In  the  horizontal  directions,  a 
unit  cell  consists  of  a  f-point  in  the  center  (fij),  a  U-point  to  its  left 

(Ui  j),  and  a  V-point  to  its  bottom  (V^  j).  In  the  vertical  direction,  the 
vertical  velocities  are  computed  at  the  "full"  giid  points.  Horizontal 
velocities,  temperature,  salinity,  and  density  are  computed  at  the  "half"  grid 
points  (half  grid  spacing  below  the  full  points). 

23.  Two  arrays,  each  of  dimension  (IMAX,  JMAX) ,  are  used  to  index  the 
grid  cells.  The  array  NS  indicates  the  condition  of  the  left  and  right  cell 
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NS 


MS 


NS 


MS 


3 

m 


Figure  4.  Computational  arrays 

boundaries,  while  the  array  MS  denotes  the  condition  of  the  top  and  bottom 
cell  boundaries  (Figure  4) . 


Numerical  Solution  Algorithm 

24.  Finite  differences  are  used  to  replace  derivatives  in  the  governing 
equations,  resulting  in  a  system  of  linear  algebraic  equations  to  be  solved  in 
both  the  external  and  internal  modes. 

25.  The  external  mode  solution  consists  of  the  surface  displacement 
and  vertically  integrated  contravariant  unit  flows  U  and  V  .  All  of  the 
terms  in  the  transformed  vertically  averaged  continuity  equation  are  treated 
implicitly  whereas  only  the  water-surface  slope  terms  in  the  transformed 
vertically  averaged  momentum  equations  are  treated  implicitly.  If  the 
external  mode  is  used  as  purely  a  vertically  avera^’^'d  model,  the  bottom 
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friction  is  also  treated  implicitly.  Those  terms  treated  implicitly  are 
weighted  between  the  new  and  old  time-steps.  The  resulting  finite  difference 
equations  are  then  factored  such  that  a  ^-sweep  followed  by  an  fy-sweep  of  the 
horizontal  grid  yields  the  solution  at  the  new  time-step. 

26.  Writing  Equations  11-13  as 


au  H  .  ar  .  „ 
St  in 


.  H  „  ar  .  „ 
-55  -  " 


the  sweep  is 


l-sweep  » ri, .  (u;.,.,  -  u;,)  = 


-  -««) 


-  ^  K.J.I  -  O 


j2  \fil  ii-l.J/ 


(l-tf) 


(rij  -  fl-x.a)  +  At  M“ 
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A  typical  value  of  0  of  0.55  yields  stable  and  accurate  solutions. 

27.  The  internal  mode  consists  of  computations  from  Equations  21-25  for 
the  three  velocity  components  u  ,  v  ,  anS  w  ,  salinity,  and  temperature. 

The  only  terms  treated  implicitly  are  the  vertical  diffusion  terms  in  all 
equations  and  the  bottom  friction  and  surface  slope  terms  in  the  momentum 
equations.  Values  of  the  water-surfac^ elevations  from  the  external  mode  are 
used  to  evaluate  the  surface  slope  terms  in  Equations  21  and  22.  As  a  result, 
the  extremely  restrictive  speed  of  a  free-surface  gravity  wave  is  removed  from 
the  stability  criteria.  Roache's  second  upwind  differencing  is  used  to  repre¬ 
sent  the  convective  terms  in  the  momentum  equations,  whereas  a  spatially 
third-order  scheme  developed  by  Leonard  (1979)  called  QUICKEST  is  used  to 
represent  the  advective  terms  in  Equations  24  and  25  for  salinity  and 
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temperature,  respectively.  For  example,  if  the  velocity  on  the  right  face  of 
a  computational  cell  is  positive  then  with  QUICKEST  the  value  of  the  salinity 
used  to  compute  the  flux  through  the  face  is 


Sr  = 


(Si, 


1 

T 


2 


1 

■? 


(41) 


The  more  interested  reader  is  referred  to  the  paper  by  Leonard  (1979)  . 

28.  It  should  be  noted  that  once  the  u  and  v  velocity  components 
are  computed  they  are  slightly  adjusted  to  ensure  the  conservation  of  mass. 
This  is  accomplished  by  forcing  the  sum  of  u  over  the  vertical  to  be  the 
vertically  averaged  velocity  U/H  and  the  sum  of  v  over  the  vertical  to 
equal  V/H  ,  where  H  is  the  total  water  depth. 

Turbulence  Parameterization 

29.  Vertical  turbulence  is  handled  by  using  the  concept  of  eddy  viscos¬ 
ity  and  diffusivity  to  represent  the  velocity  and  density  correlation  terms 
that  arise  from  a  time  averaging  of  the  governing  equations .  These  eddy  coef¬ 
ficients  are  computed  from  mean  flow  characteristics  using  a  simplified 
second-order  closure  model  originally  developed  by  Donaldson  (1973).  The 
closure  model  has  been  further  developed  and  applied  to  various  types  of  flows 
by  Lewellen  (1977)  and  Sheng  (1982,  1986).  A  discussion  of  the  implementation 
of  the  turbulence  model  taken  from  Sheng*  follows.  For  more  details  the 


*  Sheng,  Y.  P.  1990.  "A  Simplified  Second  Order  Closure  Model  of  Turbulent 
Transport,"  unpublished  paper  prepared  for  US  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 
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interested  reader  should  refer  to  these  references  and  to  Johnson  et  al. 
(1991). 

30.  Assuming  local  equilibrium  of  turbulence,  i.e.,  there  is  no  time 
evolution  or  spatial  diffusion  of  the  second-order  correlations,  an  equation 
relating  the  turbulent  kinetic  energy  and  the  macroscale  of  turbulence  to  the 
mean  flow  shear  and  stratification  (given  by  the  Richardson  number  Ri)  can  be 
derived  as 


SA^b^sQ*  +  A[(bs  +  3b  +  7b2s)Ri  -  Abs(l  -  2b)]Q2 
+  b(s  +  3  +  4bs)Ri^  +  (bs  -  A)(l  -  2b)Ri 
=  0 


where 

b  -  0.125 
s  -  1.8 
A  -  0.75 


and 


Q  = 


q 

AfJ[d\i/dz)^  +  {dv/dz)^ 


In  the  above  expression,  q  is  defined  as 


(43) 


=  (u'  u' 


V  V' 


w'w') 


1/2 


and  A  is  the  macroscale  of  turbulence .  u'  ,  v'  and  w'  are  the  turbulent 
velocity  fluctuations  and  the  overbar  indicates  time  averaging. 

31.  It  can  also  be  shown  that  the  following  relations  hold 


u'w' 


q 


w'  w' 


(44) 
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v'w'  = 


(45) 


1  + 


u 

u 


w'w' 


q^b 


(1  +  A) 
Q^d  -  «) 


+  ti) 


w'w' 


(46) 


where 


u> 


and 


u 


Ci> 

u~ 

1  -  'BS 


(47) 


(48) 


Thus,  after  the  velocity  shear  and  flow  stratification  are  determined,  q  can 


be  computed  from  Equations  42  and  43.  w'w'  is  then  determined  from 


^  -q^b 
|(1  -  2u>) 


(49) 


Finally,  after  A  is  prescribed,  u'w'  and  v'w'  can  be  computed  from  Equa¬ 
tions  44  and  45  and  the  vertical  eddy  coefficients  can  be  determined  from 


K  = 


A  +  (*> 
A(1  -  «) 


w^ 


(50) 
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(51) 


-p'w' 

dp 

-31 


A 

q 


bs 

(bs  -  w)A 


w'  w' 


In  addition  to  setting  A  -  0.65z  near  boundaries,  three  basic  constraints 
are  used  to  compute  A  at  a  vertical  position  z 


^  —  Qcut  ^g>gBax/2) 


(52) 


(53) 

(54) 


where  N  is  the  Brunt-Vaisala  frequency.  Equation  54  states  that  A  is  less 
than  a  fraction  of  the  spread  of  turbulence  as  measured  by  the  distance 
between  the  location  of  a  maximum  q^  to  where  q^  is  equal  to  25  percent  of 
the  maximum.  The  coefficient  the  order  of  0.15  to  0.25. 
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PART  III:  STRUCTURE  OF  THE  CHESAPEAKE  BAY  3D 
HYDRODYNAMIC  COMPUTER  MODEL 


32.  There  are  several  subroutines  in  the  Chesapeake  Bay  Model  plus  a 
main  program.  Subroutines  governing  model  setup  are  called  from  the  main 
program  while  subroutines  governing  the  computations  are  called  from  sub¬ 
routine  CH3DM2.  Each  subroutine  is  listed  below  with  a  description  of  its 
function.  Entry  points  in  particular  subroutines  are  also  noted. 

VCH3D  The  main  program. 

CH3D1R  Reads  data  from  FILE  4  (see  Appendix  B) .  Water  depths  at  the  cen¬ 
ter  of  each  computational  cell  are  read  from  FILE  12  (see  Appen¬ 
dix  C) .  Various  constants  are  computed,  the  layer  thicknesses  are 
set,  and  the  array  KM(I,J)  indicating  the  bottom  layer  of  each 
column  is  set. 


CH3DTR  Reads  (x,y)  coordinates  of  the  boundary-fitted  grid  from  FILE  15. 

The  coordinates  are  then  multiplied  by  the  scale  factor,  XMAP,  and 
divided  by  XREF  to  make  them  nondimens ional .  Subroutine  BJINTR  is 
called  to  provide  the  coordinate  derivatives  needed  to  compute  the 
metrics  of  the  transformation. 


BJINTR 

CH3DIH 

CH3DND 

CH3DII 

CH3DIF 

CH3DIV 


Computes  various  coordinate  derivatives  and  sets  the  water  depth 
HU(I,J)  and  HV(I,J)  on  the  faces  of  each  computational  cell. 

Prints  water  depth  arrays  if  requested  by  input  data.  Also,  the 
water  depths  are  made  nondimens ional  by  dividing  by  ZREF. 

Computes  many  constants  as  nondimens ional  quantities.  These  in¬ 
clude  the  Eckman  number,  Rosby  number,  time-step,  etc. 

Sets  up  the  arrays  that  indicate  the  types  of  boundaries  comprising 
a  computational  cell.  In  addition,  arrays  controlling  the  computa¬ 
tion  of  the  convective  terms  in  the  momenttun  equations  and  the 
averaging  of  water-surface  elevations  at  points  other  than  the  cell 
center  are  set  up. 

Initializes  various  variables  for  a  cold  run  and  creates  time  files 
for  tide,  salinity,  etc.  The  hot  start  capability  is  not 
operational. 

The  arrays  created  in  CH3DII  concerning  the  averaging  of  water- 
surface  elevations  contain  logical  values .  Those  arrays  are  used 
in  this  subroutine  to  create  arrays  containing  a  numerical  value. 
These  arrays,  i.e.  AFVI(I,J),  etc.,  are  used  to  control  the  averag¬ 
ing  of  not  only  water-surface  elevations  but  other  variables  as 
well . 
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CH3DWS  Controls  the  reading  of  either  wind  speed  or  wind  stress.  If  the 

wind  speed  is  read,  the  stress  is  computed  from  Garrett's  equation. 
ENTRY  CH3DWT  controls  the  time-varying  reads  and  computations . 

CH3DTD  Computes  a  sinusoidal  tide  from  amplitude,  phase,  and  offset  values 
input  in  CH3DIR.  ENTRY  CH3DTI  controls  the  reading  of  time-varying 
tabular  tide  data. 

33.  The  subroutines  above  are  called  from  VCH3D  in  the  sequence  given. 

Before  calling  CH3DM2 ,  which  controls  the  computations,  the  initial  salinity 

field  is  read  from  FILE  74.  The  initial  temperature  field  is  read  from 

FILE  17. 

CH3DH2  Final  subroutine  called  from  VCH3D.  All  subroutines  controlling 
the  actual  3D  computations  are  called  from  this  subroutine  in  the 
order  as  they  appear  below. 

CH3DDP  Computes  total  water  depths  and  the  thickness  of  the  top  layer  from 
the  latest  water-surface  elevation  field.  ENTRY  CH3DDM  sets  total 
water  depths  at  the  Intermediate  time  level,  and  ENTRY  CH3DDN  sets 
total  water  depths  and  the  top  layer  thickness  at  the  previous  time 
level . 

CH3DTK  Reads  equilibrium  temperatures  and  surface-heat  exchange  coeffi¬ 
cients  from  FILE  19  and  then  casts  them  into  nondimens ional  form. 
ENTRY  CH3DTB  controls  the  time-varying  read  and  interpolation. 

CH3DRI  Reads  river  inflows  from  FILE  13.  ENTRY  CH3DRV  controls  the  time- 
varying  read  and  interpolation. 

CH3DSAI  Reads  salinity  and  temperatures  at  tidal  boundaries  from  FILE  76. 

ENTRY  CH3DSAV  controls  the  time-varying  reads,  interpolation,  and 
conversion  to  nondimens ional  form. 

CH3DTEI  Reads  temperatures  at  river  inflow  boundaries  from  FILE  78.  ENTRY 
CH3DTEV  controls  the  reading  of  time— varying  temperatures  and 
interpolation. 

34.  At  this  point  the  loop  over  time  is  entered  in  CH3DM2  and  each  of 

the  subroutines  below  is  called  each  time-step. 

CH3DDE  Computes  the  water  density  using  Eckert's  equation.  The  baroclinic 
terms  in  the  momentum  equations  are  then  evaluated. 

CH3DED  Sets  up  arrays  that  are  then  used  in  SUBROUTINE  CH3DEZ  for  the 
computation  of  the  eddy  viscosity  and  eddy  diffusivity 
coefficients . 

CH3DWT  ENTRY  in  CH3DWS  for  reading  time-varying  wind  data. 

CH3DTB  ENTRY  in  CHIDTK  for  reading  time-varying  equilibrium  temperature 

and  heat-exchange  coefficients. 
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CH3DRV 

CH3DTEV 

CH3DTI 

CH3DTD 

CH3DSAV 

CH3DDN 

CH2DXY 

CH3DDP 

CH3DXYZ 

CH3DDI 

CH3DSA 

CH3DTE 

CH3D0T 

CH3DC1 

CH3DC2 

CH3DC3 

CH3DC4 


ENTRY  in  CH3DRI  for  reading  time-varying  river  flows. 

Entry  in  CH3DTEI  for  reading  time-varying  temperatures  at  river 
inflow  boundaries . 

ENTRY  in  CH3DTD  for  reading  time-varying  tabular  tide  data. 

If  a  sinusoidal  tide  is  imposed  this  subroutine  is  called  instead 
of  ENTRY  CH3DTI. 

Entry  in  CH3DSAI  for  reading  time-varying  salinity  and  temperatures 
at  tidal  boundaries. 

ENTRY  in  CH3DDP  for  assigning  total  water  depths  at  the  previous 
time  level . 

Computes  the  vertically  averaged  flow  field  from  the  vertically 
averaged  equations  of  motion. 

Using  the  water-surface  field  computed  in  CH2DXY,  computes  the 
total  water  depths . 

Computes  the  3D  velocity  field.  Mass  conservation  is  ensured  by 
forcing  the  vertical  sum  of  the  horizontal  components  of  the  3D 
velocity  to  match  the  vertically  averaged  values  computed  in 
CH2DXY. 

Computes  the  convective  and  diffusion  terms  in  the  momentum  equa¬ 
tions  using  the  most  recent  computations  from  CH3DDP  and  CH3DXYZ. 
These  terms  are  then  employed  at  the  next  time  step  in  CH2DXY  and 
CH3DXYZ. 

Computes  the  salinity  field. 

Computes  the  temperature  field. 

Controls  the  output  printed  and/or  written  to  files  for  plotting. 
Output  can  be  requested  in  terms  of  the  transformed  nondimens ional 
variables  or  the  physical  dimensional  variables.  Generally,  values 
for  the  physical  variables  in  dimensional  form  are  more  meaningful. 
If  physical  velocities  and  dimensional  output  are  requested  by 
setting  IDIM  -  1,  SUBROUTINE  CH3DC1  is  called  with  ENTRIES  CH3DC2, 
CH3DC3,  CH3DC4,  CH3DC5,  CH3DC6,  CH3DC7,  CH3DC8,  CH3DC9 ,  CH3DCA, 
CH3DCB,  CH3DCC,  CH3DCD,  CH3DCE.  Each  is  described  below. 

Provides  dimensional  water— surface  elevations. 

Provides  dimensional  physical  vertically  averaged  unit  flow  rate  in 
x-direction. 

Provides  dimensional  physical  vertically  averaged  unit  flow  rate  in 
y-direction. 

Provides  dimensional  physical  components  of  wind  shear  stress. 
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CH3DC5 

CH3DC6 

CH3DC7 

CH3DC8 

CH3DC9 

CH3DCA 

CH3DCB 

CH3DCC 

CH3DCD 

CH3DCE 

35. 

generating 
FILE  22 

FILE  23 

FILE  24 

FILE  25 

FILE  26 

FILE  31 

FILE  32 


Provides  dimensional  physical  horizontal  components  of  3D  velocity. 

Provides  dimensional  physical  vertical  component  of  3D  velocity. 

Provides  dimensional  physical  components  of  the  vertically  averaged 
velocity. 

Provides  dimensional  water  density. 

Provides  dimensional  physical  horizontal  components  of  3D  velocity 
averaged  over  some  period  of  time. 

Provides  dimensional  physical  horizontal  components  of  3D  velocity 
at  the  center  of  a  cell. 

Provides  dimensional  water  temperature. 

Provides  dimensional  physical  horizontal  components  of  3D  velocity 
at  the  center  of  a  cell  averaged  over  some  period  of  time. 

Provides  dimensional  vertical  eddy  viscosity. 

Provides  dimensional  vertical  eddy  diffusivity. 

In  SUBROUTINE  CH3D0T,  the  following  files  are  created  for  use  in 
either  time  series  plots,  vector  plots,  or  contour  plots. 

For  generating  time  series  plots  of  dimensional,  Cartesian  verti¬ 
cally  Integrated  velocities  (x  &  y  direction)  at  cell  faces  at 
specified  horizontal  locations. 

For  generating  time  series  plots  of  dimensional,  Cartesian  verti¬ 
cally  averaged  velocities  (x  &  y  direction)  at  cell  faces  at 
specified  horizontal  locations. 

For  generating  time  series  plots  of  nondimens ional  contravariant 
vertically  Integrated  velocities  at  cell  faces  at  specified 
horizontal  locations. 

For  generating  time  series  plots  of  dimensional  x  &  y-direction 
components  of  velocity  at  cell  faces  at  specified  horizontal 
locations  in  all  layers . 

For  generating  time  series  plots  of  dimensional  x  &  y-direction 
components  of  velocity  at  cell  centers  at  specified  horizontal 
locations  in  all  layers. 

For  generating  time  series  plots  of  salinity  at  specified 
horizontal  locations  in  all  layers. 

For  generating  average  salinity  plots  over  some  specified  time 
interval  at  specified  horizontal  locations  in  all  layers. 
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FILE  33 

FILE  34 

FILE  35 

FILE  36 

FILE  40 
FILE  41 

FILE  42 

FILE  43 

FILE  44 

FILE  45 


For  generating  time  series  plots  of  temperatures  at  specified 
horizontal  locations  in  all  layers. 

For  generating  time  series  plots  of  vertical  eddy  viscosity  at 
specified  horizontal  locations  in  all  layers. 

For  generating  time  series  plots  of  vertical  eddy  diffusivity  at 
specified  horizontal  locations  in  all  layers. 

For  generating  time  series  plots  of  density  at  specified  horizontal 
locations  in  all  layers. 

Geometry  of  basin  needed  to  generate  snapshot  plots. 

For  generating  snapshot  velocity  vector  plots  (nondimensional ,  con- 
travariant)  at  particular  points  in  time  in  all  layers . 

For  generating  snapshot  velocity  vector  plots  (dimensional, 
Cartesian)  at  particular  points  in  time  in  all  layers. 

For  generating  snapshot  salinity  plots  at  particular  points  in  time 
in  all  layers. 

For  generating  a  vector  plot  of  time-averaged  velocities  in  all 
layers . 

For  generating  a  time-averaged  salinity  plot  in  all  layers. 
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PART  IV:  DEMONSTRATION  OF  SETUP  OF  INPUT  FILES 


36.  To  illustrate  the  setup  of  Input  data  files  (see  Appendix  C) ,  por¬ 
tions  of  the  input  files  for  an  application  in  which  the  hydrodynamics  of  the 
Chesapeake  Bay  during  September  1983  were  simulated  are  presented.  During 
this  period  the  bay  destratified  during  a  strong  wind  event  with  a  subsequent 
restratification  as  the  wind  abated.  Rapid  cooling  of  the  surface  waters  near 
the  middle  of  the  month  contributed  to  the  destratification  process .  Results 
from  this  application  are  presented  by  Johnson,  et  al.  (1991). 

Basic  Control  Data 

37.  Input  data  controlling  the  computations  are  provided  in  FILE  4  and 
are  listed  in  Appendix  D.  These  data  are  the  response  to  the  read  statements 
presented  in  Appendix  B.  As  can  be  seen,  the  computational  time-step  is 

600  sec  with  a  total  of  3850  time-steps  simulated.  Both  temperature  and 
salinity  computations  are  made  with  the  initial  fields  frozen  for  720  time- 
steps.  As  can  be  seen  from  Figure  2,  there  are  19  river  boundaries.  The  base 
flow  for  each  is  set  to  zero;  however,  the  time-varying  river  inflow  at  the 
fall  line  is  multiplied  by  a  factor  which  varies  from  river  to  river.  This 
factor  is  the  ratio  of  the  below  fall  line  drainage  area  to  that  upstream  of 
the  fall  line.  There  are  two  tidal  boundaries.  Time-vairying  tabular  data  are 
prescribed  for  the  tides  since  ITIDE  -  -1. 

Water  Depths 

38.  The  water  depth  field  read  from  FILE  12  is  shown  in  Appendix  E. 
These  depths  are  adjusted  in  the  code  to  reflect  the  fact  that  the  thickness 
of  the  top  layer  is  initially  213.4  cm,  whereas  all  other  layers  are  152.4  cm 
thick.  These  water  depths  are  applied  at  the  center  of  each  of  the  computa¬ 
tional  cells  shown  in  Figure  2.  They  were  obtained  from  a  bathymetric  map  of 
Chesapeake  Bay  generated  by  the  Virginia  Institute  of  Marine  Sciences. 

Freshwater  Inflows 


39.  Daily  averaged  freshwater  inflows  at  the  fall  lines  of  the  James, 
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Pamunkey,  Mattaponl,  Rappahannock,  Potomac,  Patuxent,  Susquehanna,  Patapsco, 
and  Choptank  Rivers  are  given  in  Appendix  F  for  the  first  four  days  of  the 
simulation  period.  It  can  be  seen  that  after  the  day  and  hour  are  specified, 
the  location  of  the  eleven  fall  line  boundaries  and  a  discharge  for  each  are 
given.  These  data  continue  in  this  fashion  for  the  remainder  of  the  simula¬ 
tion  period.  The  inflows  at  the  lateral  locations  shown  in  Figure  2  are  set 
internally  in  the  computer  code  as  a  percentage  of  the  fall  line  flows. 

Wind  Sneed 

40.  Wind  data  are  read  from  FILE  14.  Since  IWIND  —  5,  time— varying 
wind  speed  is  read  and  the  wind  stress  is  computed  from  Equation  30.  For  the 
September  1983  application  three  wind  fields  are  read  (NWIND  -  3).  However, 
since  the  factors  WINDFX  and  WINDFY  that  multiply  the  input  winds  to  account 
for  overland  effects  are  0.0  for  the  second  wind  field,  only  the  first  and 
third  wind  fields  are  used.  These  correspond  to  data  collected  at  the  Norfolk 
International  Airport  and  the  Patuxent  River  Naval  Station.  Internal  coding 
results  in  the  wind  speed  assigned  to  each  computational  cell  being  a  linear 
interpolation  between  these  two  locations.  Cells  north  of  the  Patuxent  sta¬ 
tion  are  assigned  wind  speeds  from  the  Patuxent  station.  As  can  be  seen  in 
Appendix  G,  each  line  of  input  in  FILE  14  contains  the  day,  hour  of  the  day, 
and  the  x  and  y  components  of  the  wind  speed  from  each  of  the  three  wind 
stations.  The  wind  stresses  computed  in  the  x  and  y  components  are  then 
transformed  to  the  contravariant  components  and  by 


^  =  J  (-  ’•xY?  ^  ’■y^e) 

Hourly  wind  data  at  the  three  stations  are  presented  in  Appendix  G  for  only 
the  first  couple  of  days  in  September  1983. 

Grid  Coordinates 


41.  The  (x,y)  coordinates  of  the  computational  grid  shown  in  Figure  2 
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are  read  from  FILE  15.  File  15  was  created  from  running  the  grid  generation 
code  WESCORA.  The  form  of  these  data  for  a  few  points  Is  Illustrated  In 
Appendix  H.  The  first  line  Is  the  number  of  points  In  the  ^  and  ij  direc¬ 
tions,  respectively.  The  (x,y)  coordinates  of  each  point,  starting  with  row  1 
and  moving  from  left  to  right,  are  then  Input  as  one  point  per  line.  These 
values  correspond  to  the  arbitrary  Cartesian  system  used  to  generate  the  grid. 
They  are  multiplied  by  the  variable  XREF  to  convert  them  Into  real  world 
distances.  Points  outside  the  grid  are  assigned  a  value  of  9.0  x  10^®. 

Tabular  Tides 

42.  The  tabular  tide  data  are  read  from  FILE  16.  Each  line  of  Input 
contains  the  hour,  minutes,  and  seconds  since  the  Initiation  of  the  simulation 
along  with  the  value  of  the  water-surface  elevation  relative  to  the  top  of  the 
top  layer.  Appendix  I  Illustrates  the  form  of  these  data  In  FILE  16  for  the 
first  59  hours  of  September  1983. 

Initial  Temperature  (Salinity) 

43.  FILES  17  and  74  are  identical  in  form.  FILE  17  contains  the 
Initial  temperature  field,  whereas  FILE  74  contains  the  Initial  salinity 
field.  These  files  were  generated  by  a  computer  program  which  uses  a  few 
observed  values  to  assign  values  to  the  individual  cells  in  each  of  the  verti¬ 
cal  layers.  The  resulting  field  is  then  smoothed  in  each  horizontal  direction 
in  each  layer  before  writing  the  final  field  to  either  FILE  17  or  FILE  74. 
Appendix  J  presents  the  Initial  temperature  field  for  layer  1  and  a  portion  of 
layer  2.  Note  that  most  of  the  values  are  zero  since^the  only  active  cells  in 
layer  1  (bottom  layer)  are  in  the  deep  channel  in  the  middle  of  the  bay. 

PCflfc  £.9?  fig  Up  t  ^  ^ 

44.  Since  the  Chesapeake  Bay  model  Is  a  3D  model,  the  bottom  drag  coef¬ 
ficient  is  computed  from  Equation  32.  Thus,  FILE  18  containing  a  spatially 
varying  Manning's  n  field  is  not  needed.  In  addition,  since  the  bottom 
roughness  is  constant  FILE  20  containing  a  spatially  varying  roughness  height 
field  is  not  needed. 
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Surface  Heat  Exchange  Information 


45.  Daily  averaged  equilibrium  temperatures  and  surface  heat  exchange 
coefficients  computed  from  meteorological  data  at  the  Patuxent  Naval  Station 
are  read  from  FILE  19.  As  indicated  in  Appendix  K,  each  input  data  line  con¬ 
sists  of  the  day  and  hour  since  the  initiation  of  the  simulation  and  the  equi¬ 
librium  temperature  in  degrees  Celsius,  and  the  surface  heat  exchange 
coefficient  in  watts/mV®C. 

Tidal  Boundary  Salinity  and  Temperature 

46 .  At  tidal  boundaries ,  the  time-varying  vertical  distributions  of 
temperature  and  salinity  are  read  from  FILE  76.  As  illustrated  in  Appendix  L, 
these  data  are  grouped  in  the  following  manner.  The  day  and  hour  since  the 
initiation  of  the  simulation  is  input  on  a  line  and  followed  by  separate  lines 
of  coding  giving  the  (I,J)  location  of  the  tidal  boundary  cell  followed  by  the 
vertical  distribution  of  salinity.  A  separate  line  gives  the  same  (I,J)  loca¬ 
tion  followed  by  the  vertical  distribution  of  temperature.  These  distribu¬ 
tions  consist  of  values  in  each  active  layer  starting  with  the  surface  layer. 
Only  the  first  two  days  of  data  from  September  1983  are  presented  in 
Appendix  L  for  illustrative  purposes. 

47.  Since  ITEMP  -  1  in  FILE  4  (see  Appendix  D) ,  equilibrium  tempera¬ 
tures  are  assigned  to  the  river  inflows  and  FILE  78  does  not  exist  for  the 
September  1983  simulation. 


38 


PART  V:  SUMMARY 


48.  This  report  serves  two  purposes.  First,  a  discussion  of  the  equa¬ 
tions  solved,  turbulence  closure,  numerical  techniques  employed,  boundary 
conditions  applied,  etc.,  provides  the  interested  reader  with  theoretical 
details  concerning  the  3D  model.  Secondly,  because  sufficient  information  is 
presented  concerning  the  structure  of  the  computer  code  and  input  data 
requirements  the  report  serves  as  a  user's  guide  for  those  readers  who  may  be 
interested  in  applying  the  model.  Examples  of  the  setup  of  input  data  files 
for  the  September  1983  simulation  serve  to  strengthen  the  use  of  the  report  as 
a  user's  guide. 

49.  The  basic  equations  solved  are  first  presented  in  their  dimensional 
form  for  a  Cartesian  coordinate  system.  The  numerical  solution  of  the  govern¬ 
ing  equations  involves  an  external  mode  consisting  of  vertically  averaged 
equations  along  with  the  full  3D  internal  mode  equations.  The  Cartesian  equa¬ 
tions  for  the  two  modes  are  then  nondimensionalized  before  casting  them  into  a 
form  for  solution  on  a  boundary-fitted  coordinate  system.  In  this  transforma¬ 
tion,  only  terms  in  the  horizontal  directions  are  transformed  since  the  Carte¬ 
sian  form  is  retained  in  the  vertical  direction.  It  was  found  in  the 
Chesapeake  Bay  long-term  simulations  that  this  form  was  required  to  maintain 
the  observed  stratification  in  the  deep  channels. 

50.  A  discussion  of  boundary  conditions  reveals  that  both  the  wind 
stress  at  the  free  surface  and  the  bottom  friction  stress  are  cast  in  a  quad¬ 
ratic  form.  The  exchange  of  heat  at  the  free  surface  is  represented  using  the 
concept  of  an  equilibrium  temperature.  At  a  tidal  boundary  the  time-varying 
water  surface  and  vertical  distributions  of  salinity  and  temperature  are  re¬ 
quired,  whereas  at  river  boundaries  the  freshwater  discharge  and  the  tempera¬ 
ture  of  the  inflow  are  required. 

51.  The  method  of  finite  differences  is  used  to  solve  the  governing 
equations .  In  the  external  mode  the  vertically  averaged  continuity  equation 
is  treated  completely  implicit,  whereas  only  the  water-surface  slope  terms  are 
treated  Implicitly  in  the  momentum  equations.  The  resulting  difference  equa¬ 
tions  are  factored  into  (^,»7)  sweeps  to  yield  an  extremely  efficient  solution 
scheme.  In  the  internal  mode  equations,  the  only  terms  treated  implicitly  are 
the  vertical  diffusion  terms,  the  bottom  friction  and  the  water-surface  slope 
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terms.  The  surface  slope  terms  are  "fed"  into  the  internal  mode  computations 
from  the  external  mode . 

52 .  The  vertical  turbulence  closure  model  is  based  upon  the  assumption 
of  local  equilibrium  of  turbulence.  Thus,  all  time  derivatives  plus  the  ad- 
vection  and  diffusion  terms  in  the  Reynold's  stress  equations  are  neglected. 
This  results  in  algebraic  equations  relating  the  eddy  coefficients  to  the 
water  column  stratification,  a  turbulence  length  scale  and  the  kinetic  energy 
of  the  turbulence. 

53.  A  discussion  of  the  functions  of  all  subroutines  and  the  order  in 
which  they  are  called  provides  a  basic  understanding  of  the  flow  of  computa¬ 
tions  in  the  computer  program.  Appendices  B  and  C  provide  detailed  informa¬ 
tion  on  the  input  and  the  form  of  the  input  required  in  various  files  for 
operation  of  the  numerical  model.  A  partial  listing  of  all  the  input  files 
created  for  the  simulation  of  the  hydrodynamics  of  Chesapeake  Bay  during 
September  1983  illustrates  the  setup  of  these  input  data. 
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APPENDIX  A:  TRANSFORMED  HORIZONTAL  DIFFUSION  TERMS 


1.  When  transforming  second  derivative  terms  many  additional  terms 
arise.  For  example,  the  horizontal  diffusion  of  momentum  terms  in  Equation  14 
of  the  main  text  become: 


A1 


A3 


2.  The  transformed  horizontal  diffusion  terms  in  Equations  17  and  18  of 
the  main  text  are  somewhat  simpler  since  contravariant  velocity  components  are 
not  involved.  These  terms,  written  for  the  salinity,  are  given  below.  Iden¬ 
tical  terms  exist  for  the  temperature. 


A4 


APPENDIX  B:  LIST  OF  INPUT  DATA  IN  FILE  4 


DUMMY 

TITLE  R\m  descriptor  A80 


DUMMY 

ICHES,  IRUNYR 

ICHES  —  1  ;  For  Chesapeake  Bay  Run 

-  0  ;  Other  Areas 

IRUNYR  -  1984  Year  simulated  if 

-  1985  ICHES  —  1.  Only  used  to  control  the  generation 
1986  of  seasonal  averaged  salinity  plots 


1018 


DUMMY 

ITl,  IT2,  DT,  ISTART,  ITEST,  ITSALT  2I8,F8. 0,718 


ITl 

IT2 

DT 

ISTART  -  0 
>  0 

ITEST  -  0 
>  0 

ITSALT  - 


Starting  time  step  (  -  1  ) 

Ending  time  step 
Computational  time  step  in  Sec. 

Cold  start 

Hot  start  (not  operational) 

No  diagnostic  output 
Diagnostic  output 

Number  of  time  steps  before  salinity  and  temperature 
computations  are  initiated 


DUMMY 

IGI,  IGH,  IGM,  IGSA,  IGTE,  IGW,  IGP,  IGS,  IGC  1018 


IGI 

IGH 

IGM 

IGSA 

IGTE 

IGW 

IGP 

IGS 

IGC 


Printout  arrays  such  as  NS,  MS,  etc. 

No  printout 

Printout  all  depth  arrays 
No  printout 

Printout  Manning's  n  values 
No  printout 

Printout  initial  salinity  arrays 
No  printout 

Printout  initial  temperature  arrays 
No  printout 

Printout  the  initial  wind  shear  stress 
No  printout 

Save  cell  centered  depth  along  with  other  information  to 
FILE  40  for  snapshot  plot 
Does  not  save 

Printout  starting  S,UI,VI  for  hot  start  (not  operational) 
No  printout 

Printout  Grid  coordinates 
No  printout 


B1 


DUMMY 

IDIM,  IPA,  IPB,  ID,  JPA,  JPB,  JD,  KPA,  KPB,  KD 


1018 


IDIM 

IPA 

IPB 

ID 

JPA 

JPB 

JD 

KPA 

KPB 

KD 


DUMMY 

IPTD,  IPUI,  IPVI,  IPUA,  IPVA,  IPU3,  IPV3,  IPW3 

IPTD  ;  Time  step  interval  at  which  water  surface 
elevations  are  printed 

IPUI  ;  Time  step  interval  at  which  UI  is  printed 

IPVI  ;  Time  step  Interval  at  which  VI  is  printed 

IPUA  ;  Time  step  interval  at  which  2D  vertically  averaged  U 

is  printed 

IPVA  ;  Time  step  interval  at  which  2D  vertically  averaged  V 
is  printed 

IPU3  ;  Time  step  interval  at  which  3D  u  is  printed 

IPV3  ;  Time  step  interval  at  which  3D  v  is  printed 

IPW3  ;  Time  step  interval  at  which  3D  w  is  printed 

IPSA,  IPTE,  IPRO,  IPGA,  IPGB,  IPSS ,  IPU3A,  IPV3A 

IPSA  ;  Time  step  interval  at  which  salinity  is  printed 

IPTE  ;  Time  step  interval  at  which  temperature  is  printed 

IPRO  ;  Time  step  interval  at  which  density  is  printed 

IPGA  ;  Time  step  interval  at  which  vertical  eddy  viscosity  (GA) 
is  printed 

IPGB  ;  Time  step  interval  at  which  vertical  eddy  diffusivity  (GB) 
is  printed 

IPSS  ;  Time  step  interval  at  which  wind  shear  stress  is  printed 

IPU3A  ;  Time  step  interval  at  which  averaged  u  over  time  step 

interval  is  printed 

IPV3A  ;  Time  step  interval  at  which  averaged  v  over  time  step 
interval  is  printed 

DUMMY 

XREF,  ZREF,  UREF,  COR,  GR,  ROO,  ROR,  TO,  TR 

XREF  ;  Reference  horizontal  grid  distance 

(Maximum  horizontal  dimension  divided  by 
number  of  cells  in  that  direction  in  cm) 

ZREF  ;  Reference  depth 

(Average  typical  depth  in  cm) 

B2 


-  1  ;  Printout  in  dimensional  form 

-  0  ;  Printout  in  non-dimensional  form 

;  Beginning  I  index  for  printout 

;  Ending  I  index  for  printout 
;  I  index  printout  interval 
;  Beginning  J  index  for  printout 
;  Ending  J  index  for  printout 
;  J  index  printout  interval 
;  Beginning  K  index  (vertical)  for  printout 
:  Ending  K  index  for  printout 
;  K  index  printout  interval 


1018 


1018 


10F8 . 0 


UREF 

COR 

GR 

ROO 

ROR 

TO 

TR 


;  Reference  horizontal  velocity 
(Average  velocity  in  cm/sec) 

;  Coriolis  parameter 
;  Gravitational  acceleration  (cm/sec^) 

;  Minimum  density  expected  (gm/cc) 

;  Reference  density  (Maximum  expected)  (gm/cc) 

;  Minimum  temperature  (Celclus) 

;  Reference  temperature  (Maximum  expected)  (Celcius) 


DUMMY 

ISPAC(I),  1-1,10  1018 


ISPAC(l)  -  1  ; 

—  0  ; 
ISPAC(2)  -  1  ; 

-  0  ; 
ISPAC(3)  -  1  : 

-  0  : 
ISPAC(4)  -  1  ; 

-  0  ; 
ISPAC(5-8)  ; 
ISPAC(9)  -  0  ; 

—  2  ; 
-  3  ; 

ISPAC(IO)  ; 


Use  spatially  variable  Manning's  n  in  vertically  averaged 
model 

Use  constant  Manning's  n  in  vertically  averaged  model 
Bottom  drag  coefficient  computed  in  3D  model 
Bottom  drag  coefficient  set  to  CBF  in  3D  model 
Coriolis  ON 
Coriolis  OFF 

Variable  bottom  roughness  height 
Constant  bottom  roughness  height  given  by  BZl 
Not  used  currently 
Constant  vertical  eddy  coefficient 
Variable  vertical  eddy  coefficients 
Variable  vertical  eddy  coefficients 
computed  from  Kent-Prichard  formulation 
Not  used  currently 


DUMMY 

SSSO,  DELTAZ,  DELTAZM  10F8.0 

SSSO  ;  Initial  water  surface  elevation  relative 
to  initial  water  depth 

DELTAZ  ;  Layer  thickness  in  cm  of  all  layers  except  the  top  layer 
DELTAZM  ;  Thickness  of  top  layer  in  cm 


DUMMY 

ITB,  ZREFBN,  CTB,  BZl,  BFL,  BFCB,  MANN  I8,9F8.0 


ITB  -  1 

-  2 

ZREFBN 

CTB 

BZl 

BFL 

BFCB 

MANN 


Linear  bottom  friction  for  internal  mode 
Quadratic  bottom  friction  for  internal  mode 
Reference  height  above  bottom  in  cm 

Constant  bottom  drag  coefficient  (typical  value  of  0.003) 

Bottom  roughness  height  in  cm 

Water  depth  in  cm  below  which  the  bottom 

friction  is  increased 

Bottom  friction  coefficient  below  BFL 

Manning's  n  in  vertically  averaged  model  if  ISPAC(l)  -  0 
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DUMMY 

THETA,  THETAS 


10F8 . 0 


THETA  ;  Time  level  weighing  factor  in  external 
mode  numerical  scheme 

THETAS  ;  Time  level  weighing  factor  in  transport-diffusion  equation 


DUMMY 

ITEMP,  ISALT,  IFI,  IFD,  IFDS  1018 

ITEM?  -  2  :  Compute  temperature  (use  time-varying  temperature  as 
river  boundary  temperature) 

-  1  ;  Compute  temperature  (use  daily  equilibrium 

temperature  as  river  boundary  temperature) 

-  0  ;  No  computation  of  temperature 
ISALT  -  1  ;  Compute  salinity 

-  0  ;  No  computation  of  salinity 

IFI  -  1  ;  Compute  nonlinear  inertia  terms  in  the  momentum  equations 

-  0  ;  Do  not  compute  inertia  terms 

IFD  -  1  ;  Include  viscous  terms  in  momentum  equations 

-  0  ;  Does  not  Include 

IFDS  -  1  ;  Include  diffusion  terms  in  salinity  and  temperature 
computations 

-  0  ;  Does  not  include 


DUMMY 

SI,  S2  10F8.0 

51  ;  Constant  in  computation  of  variable  vertical 

eddy  viscosity  [GA  -  GX  (l+Sl*Ri)"“]  if  first-order 
turbulence  model  is  used 

52  ;  Constant  in  computation  of  variable  vertical 

eddy  diffusivity  [GB  -  GX  (l+S2*Ri)™2]  if  first-order 
turbulence  model  is  used 


DUMMY 

lEXP,  AVR,  AVM,  AVMl,  AHR 


lEXP  ;  Vertical  eddy  coefficient  flag 

lEXP  -  0  ;  Constant  eddy  coefficient 

—1  ;  Munk-Anderson  type  first  order  turbulence  model 
Richardson-number  dependent  eddy  coefficient 
with  length  scale  linearly  increasing  from  the 
bottom  and  surface 

—2  ;  Munk-Anderson  type  first  order  turbulence  model 
Richardson-number  dependent  eddy  coefficient 
with  length  scale  linearly  increasing  from  the 
bottom  to  the  surface 


—3 

AVR 

AVM 


Second  order  turbulence  model 

Reference  vertical  eddy  coefficient  (cm^/sec) 

Minimum  vertical  eddy  viscosity  (cmVsec) 


1I8,8F8.0 


AVMl  ;  Minimum  vertical  eddy  diffusivity  (cmVsec) 

AHR  ;  Reference  horizontal  eddy  diffusivity  or  viscosity  (cmVsec) 


DUMM\ 

FMl,  FM2,  ZTOP,  SLMIN,  QQMIN  10F8.0 

FMl  ;  Parameter  in  Richardson-number  dependent 

eddy-viscosity  (see  definition  of  SI) 

FM2  ;  Parameter  in  Richardson-number  dependent 

eddy-diffusivity  (see  definition  of  S2) 

ZTOP  ;  Distance  between  the  top  of  the  computational  domain 
and  the  free  surface.  Used  in  computing  turbulence 
length  scale  (cm) 

SUflN  ;  Minimum  value  of  turbulence  macroscale  (cm) 

QQMIN  ;  Minimum  value  of  turbulent  kinetic  energy  (gm/cm/sec^) 


DUMMY 

ICUT,  KSM,  KSL,  QCUT,  GAMAX,  GBMAX,  FZS  3I8,7F8.0 


ICUT  -  0  ;  Eddy  coefficients  constant  below  halocline 
-  1  ;  Eddy  coefficients  computed  below  halocline 
KSM  ;  Number  of  times  eddy  viscosity/diffusivity  in 

turbulence  model  are  smoothed 


KSL 

QCUT 

GAMAX 

GBMAX 

FZS 


Layer  number  above  which  smoothing  of  eddy 
viscosity/diffusivity  will  be  performed 

Coefficient  in  second-order  turbulence  model  (0.15  -  0.25) 
Maximum  value  of  eddy  viscosity  (cmVsec) 

Maximum  value  of  eddy  diffusivity  (cmVsec) 

Turbulence  scale  is  not  allowed  to  exceed  the 
product  of  FZS  and  depth 


DUMMY 

NWIND,  IWIND,  TAUX,  TAUY 


NWIND 
IWIND  -  0 

-  1 
-  2 

-  3 

-  4 

-  5 

-  6 
-  7 

TAUX 

TAUY 


Number  of  stations  where  wind  data  are  input 

Steady  and  uniform  wind  stress 

Steady  and  uniform  wind  speed 

Steady  and  space  variable  wind  stress 

Steady  and  space  variable  wind  speed 

Time  variable  and  uniform  wind  stress 

Time  variable  and  uniform  wind  speed 

Time  and  space  variable  wind  stress 

Time  and  space  variable  wind  speed 

Uniform  wind  stress  in  x-dlrection  if  IWIND  -  0 

Uniform  wind  speed  in  x-direction  if  IWIND  -  1 

Uniform  wind  stress  in  y-direction  if  IWIND  ~  0 

Uniform  wind  speed  in  y-direction  if  IWIND  -  1 


2I8,8F8.0 


B5 


If  NWIND  -  0,  use  the  following  card 


DUMMY 

If  NWIND  >  0,  use  the  following  cards 
DUMMY 

WINDFX(K),  WINDFY(K) 

WINDFX(K)  ;  Adjustment  factor  of  wind  speed  in  x  direction 
WINDFY(K)  ;  Adjustment  factor  of  wind  speed  in  y  direction 

*  Repeat  NWIND  times 


DUMMY 

XMAP 

XMAP  ;  Factor  that  scales  the  (x,y)  coordinates 
created  by  GRID  Generation  Code  to  the 
real  world 


DUMMY 

NRIVER  ;  Number  of  river  boundaries 

NRIVER  -  0  ;  No  river  boundaries 

<  0  ;  River  inflows  are  steady 
>  0  ;  Time  variable  inflows 


If  NRIVER  -  0,  use  the  following  cards 


DUMMY 

DUMMY 


If  NRIVER  >  0,  use  the  following  cards 


DUMMY 

IJRDIR(K),  IJRROW(K),  IJRSTR(K) ,  IJREND(K) 


IJRDIR(K) 


IJRROW(K) 

IJRSTR(K) 

IJREND(K) 


1  ;  River  boundary  is  on  left (west) 

2  ;  River  boundary  is  on  bottom  (south) 

3  ;  River  boundary  is  on  right  (east) 
k  ;  River  boundary  in  on  top  (north) 

;  Index  of  the  row  (J)  or  column  (I) 
of  the  river  boundary 

;  Starting  I  or  J  index  of  the  river  boundary 
;  Ending  I  or  J  index  of  the  river  boundary 


*  Repeat  NRIVER  times 
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10F8.0 


10F8.0 


1018 


1018 


DUMMY 

FACTOR(K),  BASEF(K) 


10F8.0 


FACTOR(K)  ;  Factor  by  which  time-varying  inflow  is  multiplied 

BASEF(K)  ;  Additional  inflow  added  to  time-varying  inflow  (ft^/sec) 

*  Repeat  NRIVER  times 

If  NRIVER  <  0,  use  the  following  cards 

DUMMY 

IJRDIR(K),  IJRROW(K),  IJRSTR(K) ,  IJREND(K)  1018 


IJRDIR(K)  -  1 

-  2 

-  3 

-  A 

IJRROW(K) 

IJRSTR(K) 

IJREND(K) 


;  River  boundary  is  on  left  (west) 

;  River  boundary  is  on  bottom  (south) 

;  River  boundary  is  on  right  (east) 

;  River  boundary  is  on  top  (north) 

;  Index  of  the  row  (J)  or  column  (I) 
of  the  river  boundary 

;  Starting  I  or  J  index  of  the  river  boundary 
;  Ending  I  or  J  index  of  the  river  boundary 


DUMMY 

ICELL,  JCELL,  QRIVER(K,IJ)  2I8,8F8.0 

ICELL,  JCELL  ;  Coordinate  of  a  cell(I,J)  where  QRIVER  is 
prescribed 

QRIVER(K,IJ)  :  Steady  river  inflow 
*  Repeat  NRIVER  times 


DUMMY 

NBAR  1018 

NBAR  ;  Number  of  interior  barriers 


If  NBAR  -  0,  Use  the  following  card 
DUMMY 

If  NBAR  >  0,  Use  the  following  cards 
DUMMY 

IJBDIR(K),  IJBROW(K),  IJBSTR(K) ,  IJBEND(K)  1018 

IJBDIR(K)  -  1  ;  barrier  is  horizontal  (^-direction) 

-  2  ;  barrier  is  vertical  (ij-direction) 

IJBROW(K)  ;  Index  of  row  (J)  or  column  (I)  of 

interior  barrier 
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IJBSTR<K)  ;  Starting  I  or  J  index  of  interior 
barrier 

IJBEND(K)  ;  Ending  I  or  J  index  of  interior  barrier 
*  Repeat  NBAR  times 


DUMMY 

HIDE  ;  Flags  for  tidal  boundaries 

HIDE  -  0  ;  No  tidal  boundaries 

>  0  ;  Sinusoidal  boundaries 

<  0  ;  Tabular  tidal  boxmdaries 

For  ITIDE  -  0,  use  the  following  cards 

DUMMY 

DUMMY 

DUMMY 


For  sinusoidal  tidal  boundaries,  use  the  following  cards 
DUMMY 

IJLINE,  NCG 

IJLINE  ;  Number  of  tidal  boundaries 

NCG  ;  Number  of  tidal  constituents 


DUMMY 

IJDIR(I),  IJROW(I),  IJSTR(I),  IJEND(I) 


IJDIR(I) 


IJROW(I) 

IJSTR(I) 

IJEND(I) 


1  ;  Tidal  boundary  is  on  left  (west) 

2  ;  Tidal  boundary  is  on  bottom  (south) 

3  ;  Tidal  boundary  is  on  right  (east) 

4  :  Tidal  boundary  is  on  top  (north) 

;  The  index  of  the  row  (J)  or  column  (I) 
of  the  tidal  boundary 
;  Starting  I  or  J  of  the  tidal  boundary 
:  Ending  I  or  J  of  the  tidal  boundary 


*  Repeat  IJLINE  times 


DUMMY 

AMP(I),  PHA(I),  CA(I),  TP(I) 

AMP(I)  ;  Tidal  amplitude  of  the  constituent  I  (cm) 
PHA(I)  ;  Tidal  phase  of  the  constituent  I  (degrees) 

CA  (I)  ;  Constant  tidal  amplitude  added  to  AMP(I)  (cm) 

TP  (I)  ;  Tidal  period  of  the  constituent  I  (sec) 

*  Repeat  NCG  times  for  each  IJLINE. 


1018 


1018 


1018 


10F8 . 0 
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For  tabular  tidal  boundary  conditions,  use  the  following  cards 
DUMMY 

TIDFNO,  TIDBND  10F8.0 

TIDFNO  ;  Number  of  tidal  elevation  tables  entered  as  input 

TIDBND  ;  Number  of  tidal  elevation  boundaries 


DUMMY 

TIDSTR(I),  I-l, TIDFNO  1018 

TIDSTR(I)  ;  The  entry  number  in  each  tidal  elevation 
table  corresponding  to  the  starting  time 
of  the  simulation 

DUMMY 

IJDIR(I),  IJROW(I),  IJSTR(I),  IJEND(I),  TIDTYP(I) , 

TIDFNl(I),  TIDFN2(I)  418, A8, 518 

IJDIR(I)  -  1  ;  Tidal  boundary  is  on  left  (west) 

-  2  ;  Tidal  boundary  is  on  bottom  (south) 

-  3  ;  Tidal  boundary  is  on  right  (east) 

-  4  ;  Tidal  boundary  is  on  top  (north) 

IJROW(I)  ;  Index  of  the  row  (J)  or  column  (I)  of  the 

tidal  boundary 

IJSTR(I)  ;  Starting  I  or  J  index  of  the  tidal  boundary 

IJEND(I)  :  Ending  I  or  J  index  of  the  tidal  boundary 
TIDTYP(I)  -  "CONSTANT"  ;  Constant  tidal  elevation  between 

IJSTR(I)  and  IJEND(I) 

-  "INTERP  "  ;  Linear  interpolation  of  tidal 

elevation  IJSTR(I)  and  IJEND(I) 

TIDFNl(I)  ;  The  number  of  the  tidal  elevation  table  for  CONSTANT 

or  INTERP  type  of  boundaries 

TIDFN2(I)  ;  The  number  of  the  2nd  tidal  elevation  table  used  for 

Interpolation  on  "INTERP"  t3rpe  boundary 


DUMMY 

ITBRK(I),I  -  1,  20  1018 

ITBRK(I)  ;  Time  steps  at  which  information  is  written  to 
files  for  snapshot  plots 


DUMMY 

ISBRK(I),  I-l.  10  1018 

ISBRK(l)  -  1  ;  Save  3D  dimensionless,  contravariant 
velocities  at  cell  faces 
-  0  ;  Does  not  save 
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ISBRK(2)  -  1  ;  Save  3D  cell  centered  velocities 

-  0  :  Does  not  save 
ISBRK(3)  -  1  ;  Save  salinities 

-  0  ;  Does  not  save 

ISBRK(4)  -  1  ;  Save  average  3D  velocities  over  simulation  period 
at  cell  faces 

-  0  ;  Does  not  save 

ISBRK(5)  -  1  ;  Save  average  salinities  over  simulation  period 

-  0  ;  Does  not  save 

ISBRK(6-10)  ;  Not  used  currently 


DUMMY 

NSTA,  NFREQ,  NSTART  1018 

NSTA  ;  Number  of  stations  where  information  is  saved  for 
time  series  plots  of  currents 
NFREQ  ;  Time  step  interval  for  saving  currents 
NSTART  ;  Beginning  time  step  for  saving  current  information 


DUMMY 

IST(K),  JST(K),  STATID(K)  214, A48 

IST(K),  JST(K)  ;  Coordinate  (I,J)  of  a  station  where  currents  are  saved 
STATID(K)  ;  Station  descriptor 


DUMMY 

JSPAC(I),  I  -  1,  10 

JSPAC(l)  -  1  ;  Save  dimensional  Cartesian  unit  flows  UI  and  VI  on  a  file 
-  0  ;  Does  not  save 

JSPAC(2)  -  1  ;  Save  depth-averaged  dimensional  Cartesian  U  and  V 
on  a  file 


-  0 

;  Does  not  save 

JSPAC(3)  -  1 

;  Save  dimensional  Cartesian  3D 
on  a  file 

u 

and  V  at 

cell 

faces 

-  0 

;  Does  not  save 

JSPAC(4)  -  1 

;  Save  dimensional  Cartesian  3D 

u 

and  V  at 

cell 

center 

-  0 

;  Does  not  save 

JSPAC(5)  -  1 

;  Save  nondlmenslonal  Cartesian 

UI 

and  VI 

-  0 

;  Does  not  save 

JSPAC(6-10) 

;  Not  used  currently 

DUMMY 

NSTATS,  NFREQS,  NSTRTS  1018 

NSTATS  ;  Number  of  stations  where  water  surface 

elevations  are  saved  for  time  series  plots 
NFREQS  ;  Time  step  interval  for  saving  water  surface 
elevations 

NSTRTS  ;  Beginning  time  step  for  saving  water  surface  elevations 


BIO 


DUMMY 

ISTS(K),  JSTS(K),  STATS(K) 


2I4.A48 


ISTS(K),  JSTS(K)  ;  Coordinate  (I,J)  of  a  station  where  water  surface 

elevations  are  saved 
STATS (K)  ;  Station  descriptor 

Repeat  NSTATS  Times 


DUMMY 

MSTA,  MFREQ,  MSTART 

MSTA  ;  Number  of  stations  where  salinity,  temperature, 
etc . ,  are  saved  for  time  series  plots 
MFREQ  ;  Time  step  interval  for  saving  information 
MSTART  ;  Beginning  time  step  for  saving  information 


DUMMY 

ISTSA(K),  JSTSA(K),  STATSA(K) 

ISTSA(K) ,  JSTSA(K)  ;  Coordinate  (I,J)  of  a  station  where  salinities, 

etc . ,  are  saved 

STATSA(K)  ;  Station  descriptor 

Repeat  MSTA  Times 


DUMMY 


JTPAC(I),  I  -  1, 

10 

JTPAC(l)  -  1  ; 

Save  salinities  on  file  at  MFREQ 

-  0  ; 

Does  not  save 

JTPAC(2)  -  1  ; 

Save  average  salinity  over  simulation 

period  on  file 

-  0  ; 

Does  not  save 

JTPAC(3)  -  1  ; 

Save  temperatures  on  file  at  MFREQ 

-  0  : 

Does  not  save 

JTPAC(4)  -  1  ; 

Save  eddy  viscosities  on  file  at  MFREQ 

“  0  ; 

Does  not  save 

JTPAC(5)  -  1  ; 

Save  eddy  diffusivity  on  file  at  MFREQ 

-  0  ; 

Does  not  save 

JTPAC(6)  -  1  ; 

Save  density  on  file  at  MFREQ 

-  0  ; 

Does  not  save 

JTPAC(7-10)  ; 

Not  used  currently 

DUMMY 

I,J,K,  NS(I,J,K) 

1,J  ;  Horizontal  coordinates  of  cell 
K  ;  Vertical  layer 
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2I4,A48 
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NS  ;  Indicates  type  of  left/right  boundary  of  the  cell 
(See  Figure  4,  main  text).  Allows  for  changing 
Internally  computed  values.  Peculiar  to  Chesapeake 
Bay. 


DUMMY 

I.J.R,  MS(I,J,K) 

I,J  ;  Horizontal  coordinates  of  cell 
K  ;  Vertical  layer 

MS  ;  Indicates  type  of  top/bottom  boundary  of  the  cell 
(See  Figure  4,  main  text).  Allows  for  changing 
Internally  computed  values.  Peculiar  to  Chesapeake 
Bay. 
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DUMMY 

I,J,K,  NS(I,J,K) 

I,J  ;  Horizontal  coordinates  of  cell 
K  ;  Vertical  layer 

NS  ;  Indicates  type  of  left/right  boundary  of  the  cell 
(See  Figure  4,  main  text).  Allows  for  second 
change  of  Internally  computed  values.  Peculiar 
to  Chesapeake  Bay. 
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DUMMY 

I,J,K,  MS(I,J,K) 

I,J  ;  Horizontal  coordinates  of  cell 
K  ;  Vertical  layer 

MS  ;  Indicates  type  of  top/bottom  boundary  of  the  cell 
(See  Figure  4,  main  text).  Allows  for  second 
change  of  Internally  computed  values.  Peculiar 
to  Chesapeake  Bay. 
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APPENDIX  C:  LIST  OF  DATA  INPUT  FILES 


FILE  12 

The  water  depths  at  the  center  of  each  horizontal  grid  cell  relative  to 
the  zero  of  the  dattiin,  e.g.,  NGVD,  are  read  from  FILE  12.  These  data  are 
formatted  by  16F5.0  and  are  read  by  rows  with  the  first  row  being  at  the  bot¬ 
tom  of  the  transformed  plane.  The  depths  are  input  in  units  of  feet*  and 
then  converted  to  centimeters  in  the  model. 

FILE  13 

River  inflows  are  read  from  FILE  13.  These  data  are  read  first  as  a 
time  line  (DAY  and  HOUR)  formatted  by  218.  Next  the  (I,J)  location  and  dis¬ 
charge  in  cubic  feet  per  second  for  each  cell  of  each  river  boundary  are  read 
and  formatted  by  218,F8.0. 

FILE  14 

Wind  data  are  read  from  FILE  14.  These  data  are  in  the  form  of  time 
(DAY  and  HOUR)  and  the  ^  (east-west)  and  (north-south)  components  of  the 
wind  speed  in  meters  per  second  of  each  wind  field  used.  These  data  are 
formatted  by  (215 , 6F10.0) . 

FILE  15 

The  (x,y)  coordinates  of  the  Chesapeake  Bay  grid  are  read  from  FILE  15. 
This  file  was  created  from  a  run  of  the  grid  generation  code  called  WESCORA. 
It  is  an  unformatted  file  containing  NX,  NY,  X(I,J),  Y(I,J). 

FILE  16 

Tabular  tide  data  are  read  from  FILE  16.  These  data  are  in  the  form  of 
time  (HOUR,  MINUTES,  SECONDS)  and  the  water— surface  elevations  in  centimeters 
relative  to  NGVD  for  TIDFNO  points  on  a  tidal  line.  These  data  are  formatted 
by  (3I5,8F8.0). 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  4. 
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FTLE  17 

The  initial  temperature  field  in  degrees  Celsius  is  read  from  FILE  17  by 
format  (10E12.5).  This  file  is  created  from  a  few  observed  values.  The 
resulting  field  is  then  smoothed  in  each  horizontal  direction  several  times 
before  creating  FILE  17 . 

FILE  18 

If  the  model  is  run  in  a  purely  2D  vertically  averaged  mode,  a  field  of 
Hanning's  n  values  may  be  input  by  format  (20F4.0).  The  input  values  are 
multiplied  by  0.001  to  yield  the  actual  values.  They  are  input  by  rows,  with 
the  first  row  being  at  the  bottom  of  the  transformed  plane. 

FILE  19 

Daily  average  equilibrium  temperatures  in  degrees  Celsius  and  surface 
heat  exchange  coefficients  in  units  of  watts/mV°  are  read  from  FILE  19 . 

These  data  are  in  the  form  of  time  (DAY  and  HOUR),  equilibrium  temperature, 
and  heat  exchange  coefficient.  They  are  formatted  by  (215,  FIO.O,  E12.5). 

FILE  20 

If  the  model  is  run  in  3D  mode,  a  field  of  bottom  roughness  heights  in 
centimeters  may  be  input  by  format  (20F4.0)  if  ISPAC(4)  -  1.  They  are  input 
by  rows,  with  the  first  row  being  at  the  bottom  of  the  transformed  plane. 

FILE  74 

The  initial  salinity  field  in  parts  per  thousand  is  read  from  FILE  74  by 
format  (10E12.5).  This  file  is  normally  created  from  a  few  observed  values. 
The  resulting  salinity  field  is  then  smoothed  in  each  horizontal  direction 
several  times  before  creating  FILE  74.  These  data  are  read  by  rows,  with  the 
first  row  being  at  the  bottom  of  the  transformed  plane. 

FILE  76 

Time-varying  salinity  in  parts  per  thousand  and  temperatures  in  degrees 
Celsius  at  tidal  boundaries  are  read  from  FILE  76  if  salinity  and  temperature 
are  to  be  computed.  These  data  are  in  the  form  of  time  (DAY  and  HOUR)  for¬ 
matted  by  (215).  Next  the  (I,J)  location  of  each  tidal  boundary  cell  and  a 
vertical  distribution  of  salinity,  starting  from  the  top  layer  to  the  bottom 
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layer  are  read.  These  data  are  followed  by  temperature  data  using  the  same 
format  as  for  the  salinity.  They  are  formatted  by  (215,  11F6.0). 

FILE  78 

Time-varying  temperature  at  river  flow  boundaries  are  read  from  FILE  78 
if  temperatures  are  to  be  computed  and  equilibrium  temperatures  are  not  used 
as  river  flow  boundary  temperatures.  These  data  are  in  the  form  of  a  time 
(DAY  and  HO!IR)  formatted  by  (215).  Next  read  (I,J)  location  of  river  flow 
boundary  cells  and  corresponding  temperatures  starting  from  top  layer  to  bot¬ 
tom  layer.  These  data  are  formatted  by  (2I5.11F6.0) . 
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APPENDIX  D:  INPUT  DATA  IN  FILE  4  FOR 
SEPTEMBER  1983  APPLICATION 

CHESAPEAKE  BAY  ;  SEP83-RONC26.7  DAYS),  9  MAJOR  RIVERS  WITH  10  LATERAL  FLOWS, 
WIND  STRESS,  2  TABULAR  TIDES,  SALINITY  AND  TEMPERATURE 

ICHES  IRUNYR 
1  1983 


ITl 

IT2 

DT 

ISTART 

ITEST 

ITSALT 

1 

3850 

600.0 

0 

0 

720 

IGI 

IGH 

IGM 

IGSA 

IGTE 

IGW 

IGP 

IGS 

IGC 

1 

0 

0 

0 

0 

0 

0 

0 

0 

IDIM 

IPA 

IPB 

ID 

JPA 

JPB 

JD 

KPA 

KPB 

KD 

1 

1 

67 

1 

1 

34 

1 

1 

15 

2 

IPTD 

IPUI 

IPVI 

IPUA 

IPVA 

IPU3 

IPV3 

IPW3 

1440 

0 

0 

0 

0 

1440 

1440 

1440 

IPSA 

IPTE 

IPRO 

IPGA 

IPGB 

IPSS 

IPU3A 

1PV3A 

1440 

1440 

0 

0 

0 

0 

3850 

3850 

XREF 

ZREF 

UREF 

COR 

GR 

ROO 

ROR 

TO 

TR 

452628. 

1000. 

10.0 

.0001 

981.0 

1.0 

1.021 

1. 

30. 

ISPAC(l) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

1 

1 

1 

0 

2 

SSSO 

DELTAZ 

DELTAZM 

0.0 

152.4 

213.4 

ITB 

ZREFBN 

CTB 

BZl 

BFL 

BFCB 

MANN 

2 

5.0 

0.0025 

0.005 

100. 

0.008 

0.025 

THETA 

THETAS 

1.0 

1.0 

ITEMP 

I  SALT 

IFI 

IFD 

IFDS 

1 

1 

1 

0 

0 

SI 

S2 

10. 

3.33 

lEXP 

AVR 

AVM 

AVMl 

AHR 

-3 

10. 

10.0 

0.005 

100000. 

FMl 

FM2 

ZTOP 

SLMIN 

QQMIN 

-.5 

-1.5 

0. 

1. 

0.01 

ICUT 

KSM 

KSL 

QCUT 

GAMAX 

GBMAX 

FZS 

0 

1 

1 

0.15 

5000. 

5000. 

0.2 

NVIND 

IWIND 

TAUX 

TAUY 

3 

5 

0.0 

0.0 

WNDFX  WNDFY 

2.03  1.40 

0.00  0.00 

1.53  1.85 

XMAP 
570738.1 
NRIVER 
19 
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IJRDIR 

IJRROW 

IJRSTR 

IJREND 

4 

34 

2 

2 

1 

2 

31 

31 

3 

3 

24 

24 

4 

27 

6 

6 

1 

6 

22 

22 

4 

27 

8 

8 

3 

8 

25 

25 

4 

31 

17 

17 

1 

16 

19 

19 

3 

17 

26 

26 

4 

34 

22 

22 

1 

21 

26 

26 

3 

26 

21 

21 

3 

40 

15 

15 

2 

15 

34 

34 

4 

15 

38 

38 

4 

14 

62 

62 

4 

16 

48 

50 

2 

1 

37 

37 

FACTOR 

6ASEF 

1.507 

0. 
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0. 
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0. 
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0. 

1.699 

0. 

1.699 

0. 

1.699 

0. 
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0. 

1.808 

0. 
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0. 
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0. 
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0. 
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0. 
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0. 
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0. 
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0. 

1.015 

0. 

1.346 

0. 

6.602 

0. 

(  ONE  CARD  FOR  EACH  RIVER  ) 

JAMES  R. 

Appomattox  R.  at  Matoaca 
James  R-  Lat.  2 
YORK  (PAMUNKEY  R.) 

Pamunkey  R.  Lat.  1 
YORK  (MATTAPONI  R.) 

Hattaponi  R.  Lat.  1 
RAPPAHANNOCK  R. 

Rappahannock  R.  Lat.  1 
Rappahannock  R.  Lat.  2 
POTOMAC  R. 

Potomac  R.  Lat.  1 
Potomac  R.  Lat.  2 
PATUXENT  R. 

Patuxent  R.  Lat.  1 
Patuxent  R.  Lat.  2 
SUSQUEHANNA  R. 

BALTIMORE  HARBOUR ( PAT APSCO  R. ) 
CHOPTANK  R. 
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NBAR 

9 

C 

IJRDIR 

IJRROW  IJRSTR 

IJREND 

(  ONE 

CARD  FOR 

EACH  BAR 

) 

1 

9  30 

33 

2 

34  8 

8 
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-1 

TIDFNO 
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1 

2 
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1 

(2)  (3) 

(A) 
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X 
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IJEND 

TIDTYP 
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1 

5  4 
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1 

1 

2 

4  1 
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1 

1 

ITBRK(l) 
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0  0 
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0  0 
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0 

0 

0 
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(4) 

(5) 
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(7) 
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(9) 

(10) 

1 

0  0 

0 

1 

NSTA 

NFREQ  NSTART 

(CURRENT  STATIONS) 

4 

3  1 

1ST  JST 

STATID(K)  (214, A48) 

(  ONE 
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EACH  STATION  ) 

3  5 

BAY  ENTRANCE 

13  9 

WOLF  TRAP 

31  11 

MID  BAY 

43  10 
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7 

3  1 
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EACH  STATION  ) 

3  5  CHESAPEAKE  BAY  TUNNEL,  VA  (6008863) 

1  13  HAMPTON  ROADS,  VA  (6008610) 

23  21  COLONIAL  BEACH.  VA  (6005150) 
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JTPAC(l)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 

10  10  0  0 
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15 

IS 

13 

13 

15 

20 

30 

20 

10 

10 

10 

10 

0 

0 

0 

0 

0 

0 

0 

20 

0 

0 

0  0  0 


El 


13  50  20  0  0  15  35  33  23  25  15  0  0  0  15  20 

30  43  0  0  0  0  20  30  SO  20  0  0  0  25  40  20 

000000000000000  20 
50  20  0  10  000000000  20  00 

0  0  0 

13  50  15  0  0  10  30  25  20  20  15  0  0  0  0  15 

30  45  0  0  0  0  20  SO  SO  20  0  0  0  25  40  20 

000000000000000  20 
30  20  0  10  0000000  00  20  00 

0  0  0 

10  SO  IS  0  0  10  30  25  20  20  15  0  0  0  0  10 

30  45  0  0  0  0  20  50  SO  20  0  0  0  25  40  30 

30  25  20  IS  10  10  10  10  0  0  0  0  0  0  0  20 

SO  20  0  10  000000000000 

0  0  0 

10  SO  IS  0  0  10  30  23  20  15  IS  0  0  0  0  10 

25  45  0  0  0  0  20  50  SO  20  0  0  0  IS  15  IS 

IS  IS  15  15  0  0  0  0  0  0  0  0  0  0  0  20 

SO  20  00000000000000 

0  0  0 

10  20  10  0  0  10  30  25  0  IS  IS  0  0  0  0  10 

20  0  0  0  0  0  13  40  40  15  0  0  0  0  0  0 

OOOOOOOOOOOOOOOO 
0000000000000000 
0  0  0 

10  20  10  0  0  10  30  20  0  IS  IS  0  0  0  0  10 

13  0  0  0  0  0  IS  40  40  IS  0  0  0  0  0  0 

OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
0  0  0 

10  20  10  0  0  10  so  IS  0  IS  IS  0  0  0  0  10 

10  0  0  0  0  0  IS  40  40  IS  0  0  0  0  0  0 

0000000000000  000 

OOOOOOOOOOOOOOOO 
0  0  0 

10  20  10  0  0  10  so  10  0  0  0  0  0  0  0  10 

10  0  0  0  0  0  IS  IS  30  30  0  0  0  0  0  0 

OOOOOOOOOOOOO  000 

0000000000300000 
0  0  0 

0  IS  10  0  0  10  30  10  0  0  0  0  0  0  0  10 

10  0  0  0  0  0  13  13  30  30  0  0  0  0  0  0 

OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOO  000 

0  0  0 

0  IS  10  0  0  10  23  10  0  0  0  0  0  0  0  10 

10  0  0  0  0  0  13  15  15  55  0  0  0  0  0  0 

OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
0  0  0 

0  15  10  0  0  10  20  10  0  0  0  0  0  0  0  10 

10  0000000  15  55  000  000 

OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
0  0  0 

0  15  10  0  0  10  20  10  0  0  0  0  0  0  0  0 

10  0  0  0  13  IS  IS  IS  IS  55  0  0  0  0  0  0 

OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
0  0  0 
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APPENDIX  F:  RIVER  INFLOWS  IN  FILE  13 
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APPENDIX  G:  WIND  SPEED  DATA  IN  FILE  14 


Patuxent 


Norfolk  River 

Int'l  Naval  Baltimore-Washington 


Station 


Airport 


DAY 

HOUR 

E-W 

N-S 

E-W 

N-S 

E-W 

N-S 

0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0 

X 

-0.9927 

•1.1031 

-3.1208 

-1.8018 

-2.0992 

-1.3888 

0 

2 

-0.4037 

-0.1701 

-3.1208 

•1.8018 

-2.0992 

-3.3888 

0 

3 

0.0000 

0.0000 

•3.0033 

0.0000 

-3.0418 

-0.9384 

0 

k 

-0.2902 

-1.3209 

-3.0039 

0.0000 

-3.0418 

-0.3384 

0 

S 

-X.3379 

-0.7722 

-3.0039 

0.0000 

-3.0418 

-0.3384 

0 

6 

-2.2291 

-1.2070 

-2.0472 

3.1348 

-2.9348 

0.4470 

0 

7 

-3.1200 

-1.0010 

-2.0472 

3.1348 

-2.3348 

0.4470 

0 

8 

-3.3000 

-2.0392 

-2.0472 

3.1348 

-2.9348 

0.4470 

0 

9 

-3.0410 

-0.9304 

-2.7009 

2.3183 

-1.9894 

2.3881 

0 

XO 

-4.1103 

0.0000 

-2.7003 

2.3183 

-1.9834 

2.3881 

0 

IX 

-4.0390 

0.7131 

-2.7003 

2.3183 

•1.9834 

2.3881 

0 

X2 

-4.1103 

0.0000 

0.0000 

2.0992 

-0.8237 

3.3488 

0 

13 

-3.3090 

3.9439 

0.0000 

2.0392 

-0.8237 

3.9488 

0 

14 

-2.4107 

0.0003 

0.0000 

2.0392 

-0.8237 

3.3488 

0 

IS 

•2.0279 

0.3370 

0.0000 

0.0000 

-0.3978 

2.0279 

0 

16 

-2.3740 

0.0000 

0.0000 

0.0000 

-0.3978 

2.0279 

0 

17 

-2.4107 

0.0003 

0.0000 

0.0000 

-0.3978 

2.0279 

0 

18 

-2.3001 

1.9034 

0.0000 

•2.0992 

0.9927 

1.1831 

0 

19 

0.0000 

0.0000 

0.0000 

-2.0992 

0.9927 

1.1831 

0 

20 

-0.2002 

1.3209 

0.0000 

•2.0992 

0.9927 

1.1831 

0 

21 

2.0749 

1.3444 

-0.3978 

-2.0279 

0.0000 

2.0392 

0 

22 

1.9710 

1.0343 

•0.3378 

•2.0279 

0.0000 

2.0392 

0 

23 

1.2070 

2.2291 

-0.3378 

•2.0279 

0.0000 

2.0992 

0 

24 

1.0349 

1.9718 

0.0000 

0.0000 

-0.3378 

2.0279 

1 

1 

0.0010 

0.7007 

0.0000 

0.0000 

•0.3978 

2.0279 

1 

2 

0.0000 

0.0000 

0.0000 

0.0000 

-0.3978 

2.0279 

1 

3 

0.0000 

0.0000 

0.0000 

0.0000 

•0.8049 

4.3827 

1 

4 

0.0000 

0.0000 

0.0000 

0.0000 

-0.8049 

4.3827 

1 

i 

0.0000 

0.0000 

0.0000 

0.0000 

-0.8049 

4.9827 

1 

6 

0.0000 

0.0000 

0.9384 

3.0418 

o.oooo 

6.1779 

1 

7 

0.0000 

0.0000 

0.3384 

3.0418 

o.oooo 

6.1779 

1 

8 

-1.3230 

1.3774 

0.9384 

3.0418 

0.0000 

6.1779 

1 

9 

-0.4470 

2.3340 

1.9894 

2.3881 

-1.3848 

4.3337 

1 

10 

0.0000 

3.0039 

1.9894 

2.3881 

-1.3848 

4.3937 

1 

11 

0.0297 

3.3400 

1.9894 

2.3881 

•1.3848 

4.3337 

1 

12 

-0.9904 

3.0410 

2.7809 

2.3183 

-1.4089 

1.8700 

1 

13 

-0.3304 

3.0418 

2.7803 

2.3183 

-1.4089 

3.8700 

1 

14 

1.4009 

3.8700 

2.7809 

2.3163 

-1.4089 

3.8700 

1 

13 

0.0003 

2.4187 

1.2870 

2.2291 

0.0000 

4.1183 

1 

16 

0.0000 

2.3740 

1.2870 

2.2291 

0.0000 

4.1183 

1 

17 

0.4470 

2.3348 

1.2870 

2.2291 

0.0000 

4.1183 

1 

18 

0.0000 

2.0992 

0.0000 

0.0000 

-1.2323 

1.3882 

1 

19 

-0.7043 

1.9330 

0.0000 

0.0000 

-1.2323 

3.3882 

1 

20 

-0.3370 

2.0279 

0.0000 

0.0000 

-1.2323 

1.3882 

1 

21 

-0.0003 

2-.  4107 

0.0000 

0.0000 

-2.3740 

4.4382 

1 

22 

-0.3202 

1.4312 

0.0000 

0.0000 

-2.3740 

4.4982 

1 

23 

-0.3202 

1.4912 

0.0000 

0.0000 

-2.3740 

4.4982 

1 

24 

-0.9202 

1.4312 

-2.4187 

-0.8803 

-1.1348 

2.8472 

2 

X 

•0.7007 

0.0018 

-2.4187 

-0.8803 

-1.1348 

2.8472 

2 

2 

0.0000 

0.0000 

-2.4187 

-0.8803 

-1.1348 

2.8472 

2 

3 

•1.3373 

0.7722 

0.7043 

-1.9390 

-4.8374 

1.7807 

2 

4 

-1.3373 

0.7722 

0.7043 

-1.9390 

-4.8374 

1.7807 

2 

5 

•1.9710 

1.0943 

0.7043 

-1.9390 

•4.8374 

1.7807 

2 

6 

-0.0010 

0.7807 

-1.0984 

2.9023 

•1.1348 

2.8472 

2 

7 

•1.7033 

1.0290 

-1.0984 

2.9023 

-1.1348 

2.8472 

2 

8 

-0.9075 

0.3321 

-1.0384 

2.9023 

-1.1348 

2.6472 

2 

9 

-1.7033 

1.0290 

-2.9023 

-1.0364 

-1.1208 

1.8018 

2 

10 

-1.7033 

1.0290 

-2.9029 

-1.0364 

-1.1208 

1.8018 

2 

XX 

-1.7033 

1.0290 

-2.9023 

-1.0364 

-1.1298 

1.8018 

G1 


APPENDIX  H;  COORDINATES  OF  THE  COMPUTATIONAL  GRID 

IN  FILE  15 


68 

9.000C302E-t'18 

9.000000224-18 

9.00000022-418 

9.00000022-418 

9.00000022-418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

17.80000 
19.00000 

19.80000 
21.20000 
23.20000 
23.60000 
26.40000 
24.90000 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

37.63000 

37.72000 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 


33 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

3.200000 

3.000000 

2.000000 

2.000000 

4.300000 

4.400000 

4.700000 

4.000000 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

6.100000 

6.130000 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418' 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 

9.00000022418 


HI 


APPENDIX  I:  TIDE  DATA  AT  THE  OCEAN  BOUNDARY  IN  FILE  16 


Hr  min  sec 

0  0  0 

10  0 
2  0  0 

3  0  0 

4  0  0 

5  0  0 

6  0  0 

7  0  0 

8  0  0 

9  0  0 

10  0  0 

11  0  0 

12  0  0 

13  0  0 

14  0  0 

13  0  0 

16  0  0 

17  0  0 

18  0  0 

19  0  0 

20  0  0 

21  0  0 

22  0  0 

23  0  0 

24  0  0 

23  0  0 

26  0  0 

27  0  0 

26  0  0 

29  0  0 

30  0  0 

31  0  0 

32  0  0 

33  0  0 

34  0  0 

33  0  0 

36  0  0 

37  0  0 

38  0  0 

39  0  0 

40  0  0 

41  0  0 

42  0  0 

43  0  0 

44  0  0 

43  0  0 

46  0  0 

47  0  0 

48  0  0 

49  0  0 

50  0  0 

31  0  0 

52  0  0 

53  0  0 

54  0  0 

55  0  0 

56  0  0 

57  0  0 

58  0  0 

59  0  0 


Water 
Surface 
El  ■  cm 


0.00 

20.66 

42.37 

43.42 
43.89 

37.49 

20.42 
1.32 

-3.79 

-7.01 

-6.71 

11.38 

27.43 
43.42 
37.91 
63.33 

66.14 

36.69 
43.42 
29.87 
12.80 
-2.13 
-7.92 
-4.27 

3.79 

19.20 

34.44 
43.11 
47.24 

44.50 

35.97 

21.64 
7.01 

-4.27 

-7.92 

-3.49 

6.40 

23.77 

42.98 
57.61 
67.36 
67.67 

63.70 

53.64 

34.14 
13.72 
-1.52 
-4.88 
-1.32 

8.23 

24.38 

39.32 

49.38 

58.32 
55.47 
44.30 
29.57 
12.19 
-1.83 

-10.06 


II 


0 .  OOOOOE-fOO 
O.OOOOOE-fOO 
O.OOOOOE^OO 
0 .  OOOOOE-fOO 
O.OOOOOE-fOO 
0 .  OOOOOE't’OO 
0 .  OOOOOE-fOO 
O.OOOOOE^OO 
O.OOOOOE-fOO 
0 . OOOOOE+OO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
0 .  OOOOOE'fOO 
O.OOOOOE-i-00 
O.OOOOOE't-00 
O.OOOOOE-fOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOE-fOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 


APPENDIX  J:  INITIAL  TEMPERATURE  FIELD  IN  FILE  17 


O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 


O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
0 . OOOOOEf 00 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 


O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
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O.OOOOOE-fOO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
0 .  OOOOOE't’OO 
O.OOOOOE-fOO 
O.OOOOOE+00 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
O.OOOOOE-fOO 
O.OOOOOEfOO 
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O.OOOOOE-fOO 
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O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
0.24766Ef02 
0.24826Ef02 
O.OOOOOEfOO 
O.OOOOOEfOO 
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O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
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O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 


O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

0.2476SEf02 

0.24830Ef02 

O.OOOOOEfOO 
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O.OOOOOEfOO 
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O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

0.24804Ef02 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

0.24830Ef02 

0.24765Ef02 

0.24779Ef02 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 
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O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
0.24817Ef02 
0.2476SEf02 
0.24770Ef02 
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O.OOOOOEfOO 
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O.OOOOOEfOO 
O.OOOOOEfOO 
O.OOOOOEfOO 
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O.OOOOOEfOO 
O.OOOOOEfOO 
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O.OOOOOEfOO 
O.OOOOOEfOO 


O.OOOOOEfOO 

O.OOOOOEfOO 
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O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 

O.OOOOOEfOO 
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APPENDIX  K:  EQUILIBRIUM  TEMPERATURE  AND  SURFACE 
HEAT  EXCHANGE  COEFFICIENT  IN  FILE  19 


Temp 


Day 

Sr 

»c 

0 

0 

24.20000 

1 

0 

23.60000 

2 

0 

27.20000 

3 

0 

26.00000 

4 

0 

28.70000 

5 

0 

28.80000 

6 

0 

29.40000 

7 

0 

23.40000 

8 

0 

26.90000 

9 

0 

26.90000 

10 

0 

25.80000 

11 

0 

25.40000 

12 

0 

21.30000 

13 

0 

X7.X0000 

14 

0 

X8.X0000 

15 

0 

19.40000 

16 

0 

22.30000 

17 

0 

22.60000 

18 

0 

23.40000 

19 

0 

23.40000 

20 

0 

20.80000 

2X 

0 

12.90000 

22 

0 

13.40000 

23 

0 

13.70000 

24 

0 

20.80000 

29 

0 

18.20000 

26 

0 

21.30000 

27 

0 

17.30000 

28 

0 

14.30000 

29 

0 

16.90000 

30 

0 

20.80000 

31 

0 

19.90000 

Surface 

Heat 

Exchange 

Coef.  watts /m^/°C 

0.99663E-03 

0.7S046E-03 

0.47a00E>03 

0.10779E-02 

0.83172E-03 

0.9e946E-03 

0.66442E-03 

0.6MS1E-03 

0.S9989E-03 

0.70S0SE>03 

0.10086E-02 

0.102S3E-02 

0.10373E-02 

0.11042E-02 

0.9X039E>03 

0.71461E-03 

0.78392B-93 

0.70744E.03 

0.U783E-02 

0.14483E-02 

0.1S0S7E.02 

0.1U14E-02 

0.S0X90E-03 

0.33297E-03 

0.3228SE-03 

0.3X787E.03 

0.44893E-03 

Q.98936E-03 

C.X3X2XE-02 

0.X3958E-02 

0.33336E-03 

0.461Z7E-03 
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APPENDIX  L:  TIME-VARYING  VERTICAL  DISTRIBUTIONS  OF  SALINITY 
AND  TEMPERATURE  AT  THE  OCEAN  BOUNDARY  IN  FILE  76 


0  0 

S  4  28.30  28.30 

S  4  24.90  24.90 

1  4  28.30  28.30 

1  4  2S.00  25.00 

2  4  28.30  28.30 

2  4  24.90  24.90 

3  4  28.30  28.30 

3  4  24.90  24.90 

4  4  28.30  28.30 

4  4  24.90  24.90 

5  4  28.30  28.30 

5  4  24.90  24.90 

0  12 

5  4  27.60  27.60 

5  4  24.90  24.90 

1  4  27.60  27.60 

1  4  24.90  24.90 

2  4  27.60  27.60 

2  4  24.90  24.90 

3  4  27.60  27.60 

3  4  24.90  24.90 

4  4  27.60  27.60 

4  4  24.90  24.90 

5  4  27.60  27.60 

5  4  24.90  24.90 

1  0 

5  4  28.20  28.20 

5  4  23.20  25.20 

1  4  28.20  28.20 

1  4  25.20  23.20 

2  4  28.20  28.20 

2  4  25.20  25.20 

3  4  28.20  28.20 

3  4  25.20  25.20 

4  4  28.20  28.20 

4  4  23.20  25.20 

3  4  28.20  28.20 

5  4  23.20  23.20 

1  12 

3  4  28.10  28.10 

3  4  23.00  23.00 

1  4  28.10  28.10 

1  4  23.00  23.00 

2  4  28.10  28.10 

2  4  23.00  25.00 

3  4  28.10  28.10 

3  4  23.00  23.00 

4  4  28.10  28.10 

4  4  23.00  23.00 

5  4  28.10  28.10 

5  4  23.00  23.00 

2  0 

3  4  28.10  28.10 

5  4  24.80  24.80 

1  4  28.10  28.10 

1  4  24.80  24.80 

2  4  28.10  28.10 

2  4  24.80  24.80 

3  4  28.10  28.10 


28.30  28.30 
24:90  24.90 
28.30  28.30  28.30 
23.00  25.00  25.00 
28.30  28.30  28.30 
24.90  24.90  24.90 
28.30  28.30  28.30 
24.90  24.90  24.90 
28.30  28.30 
24.90  24.90 
28.30  28.30 
24.90  24.90 

27.60  27.60 
24.90  24.90 
27.60  27.60  27.60 
24.90  24.90  24.90 
27.60  27.60  27.60 
24.90  24.90  24.90 
27.60  27.60  27.60 
24.90  24.90  24.90 
27.60  27.60 
24.90  24.90 
27.60  27.60 
24.90  24.90 

28.20  28.20 

23.20  23.20 

28.20  28.20  28.20 

23.20  23.20  23.20 

28.20  28.20  28.20 

23.20  23.20  23.20 

28.20  28.20  28.20 

23.20  23.20  23.20 

28.20  28.20 

23.20  23.20 

28.20  28.20 
23.20  23.20 

28.10  28.10 
23.00  23.00 
28.10  28.10  28.10 
23.00  23.00  23.00 
28.10  28.10  28.10 
25.00  23.00  25.00 
28.10  28.10  28.10 
23.00  23.00  23.00 
28.10  28.10 
23.00  23.00 
28.10  28.10 
23.00  25.00 

28.10  28.10 
24.80  24.80 
28.10  28.10  28.10 
24.80  24.80  24.80 
28.10  28.10  28.10 
24.80  24.80  24.80 
28.10  28.10  28.10 


29.30  29.50  29.50 
23.00  23.00  23.00 
29.50  29.50  29.50 
23.00  23.00  23.00 
29.30  29.30 
23.00  23.00 


28.80  28.80  28.80 
22.90  22.90  22.90 
28.80  28.80  28.80 
22.90  22.90  22.90 
28.80  28.80 
22.90  22.90 


29.70  29.70  29.70 
23.80  23.80  23.80 
29.70  29.70  29.70 
0.80  23.80  23.80 
29.70  29.70 
23.80  23.80 


29.80  29.80  29.80 
23.00  23.00  23.00 
29.80  29.80  29.80 
23.00  23.00  23.00 
29.80  29.80 
23.00  23.00 


29.40  29.40  29.40 

24.40  24.40  24.40 

29.40  29.40  29.40 

24.40  24.40  24.40 

29.40  29.40 


29.50  29.50  29.50 
23.00  23.00  23.00 
29.50  29.50  29.50 
23.00  23.00  23.00 


28.80  28.80  28.80 
22.90  22.90  22.90 
28.80  28.80  28.80 
22.90  22.90  22.90 


29.70  29.70  29.70 
23.80  23.80  23.80 
29.70  29.70  29.70 
23.80  23.80  23.80 


29.80  29.80  29.80 
23.00  23.00  23.00 
29.80  29.80  29.80 
23.00  23.00  23.00 


29.40  29.40  29.40 

24.40  24.40  24.40 

29.40  29.40  29.40 

24.40  24.40  24.40 


LI 


